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Abstract. We derive closed formulas for mean values of all powers of r in nonrelativistic and 
1^ , relativistic Coulomb problems in terms of the Hahn and Chebyshev polynomials of a discrete variable. 

' A short review on special functions and solution of the Coulomb problems in quantum mechanics is 

^ , given. 

(N 

Ph' 1. Introduction 



A basic problem in quantum theory of the atom is the problem of finding solutions of the non- 



^ ■ relativistic Schrodinger and relativistic Dirac wave equations for the motion of electron in a central 



attractive force field. The only atom for which these equations can be solved explicitly is the sim- 
plest hydrogen atom, or, in general, the one electron hydrogen-like ionized atom with the charge of 
the nucleus Ze; this is a classical problem in quantum mechanics which is studied in great detail; 
r--- [ see, for example, [1], [17], [16], [29], [39], [U], [51], [51], [61] and references therein. Compari- 
^ son of the results of theoretical calculations with experimental data provides accurate tests of the 
^ validity of the quantum electrodynamics [T^, [H], [82]. Explicit analytical solutions for hydrogen- 
', like atoms can be useful as the starting point in approximate calculations of more sophisticated 
O ' quantum-mechanical systems. 

i> : 

O ■ The Schrodinger equation for the hydrogen atom can be solved explicitly in the spherical and 

K*" ' parabolic system of coordinates [T7], [51] and it can be shown that the Hahn polynomials connect 

^ ! the wave functions derived in the spherical and parabolic coordinates; see [70] and references therein 

^ \ for more details. In the present paper we discuss another connection with the classical polynomials 

■ - - ' — we derive closed formulas for the mean values 

(rP) = / dv (1.1) 



R3 

of all admissible powers of r for both nonrelativistic and relativistic Coulomb wave functions in the 
spherical coordinates in terms of the Hahn and Chebyshev polynomials of a discrete variable; see, 
for example, [2], [3], [1], [5], [6], [7], [8], [36], [12], [M], [19], [55], [56], [57], [71], and [75], [76], [77] 
for an introduction to the theory of the classical orthogonal polynomials. Special cases p = 0, 1, 2 
of fll.ip give the normalization of the wave functions, the average distance between the electron and 
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the nucleus, and the mean square deviation of the nucleus-electron separation, respectively. Special 
cases p = —1, —2, —3 of these matrix elements are important in calculations of the energy levels 
by the virial theorem, the fine structure of the energy levels in Pauli's theory of the spin, radiative 
corrections and Lamb shift in hydrogen-like atoms; see [1], [IT], [16], [17], [51], [M], ^82] for more 
details. The general formulas for the corresponding relativistic Coulomb matrix elements may also 
be important in developing the theory of spectra of heavy ions for large values of Z on the basis of 
the methods of quantum field theory [19], [20], [78] . 

It is worth noting that exact solutions of certain time-dependent Schrodinger equations are found 
in [52] and [53] . An interesting area of research in physics in general is a problem of discretization of 
the space-time continuum and the concept of the fundamental length [16], [60] and, in particular, 
the discretization of the classical Maxwell, Schrodinger and Dirac equations; see [37] and [72] for 
some solutions of the discrete wave. Maxwell and Dirac equations. 

The paper is organized as follows. In the next section, among other things, we evaluate an integral 
of the product of two Laguerre polynomials in terms of the Hahn polynomials, which gives a "master 
formula" for evaluation of the matrix elements (11.11) for the nonrelativistic and relativistic hydrogen- 
like atoms in sections 3 and 4, respectively. Some special cases are given explicitly and evaluation 
of the effective electrostatic potential in the hydrogen-like atoms is discussed as an application. 
Sections 5 and 6 are written in order to make our presentation as self-contained as possible — they 
contain a short review of Nikiforov and Uvarov's approach to the theory of special functions of 
mathematical physics and a detailed solution of the wave equation of Dirac for Coulomb potential, 
respectively. We follow [57] with somewhat different details; for example, in Section 5 we give a 
different proof of the main integral representation for the special functions of hypergeometric type 
and discuss the power series method; in Section 6 we construct the spinor spherical harmonics and 
separate the variables in the spherical coordinates in detail before solving the radial equations. In 
Section 7 we discuss a more general version of the method of separation of the variables for Dirac- 
type systems. The relativistic Coulomb wave functions are not well known for a "general audience" 
and this discussion might be useful for the reader who is not an expert in theoretical physics; our 
paper is written for those who study quantum mechanics and would like to see more details than in 
the classical textbooks pp, [T7j, [IH], [H]; it is motivated by a course in quantum mechanics which 
one of the authors (SKS) has been teaching at Arizona State University for many years. Appendix 
contains some formulas which are widely used throughout the paper. 

We use the absolute cgs system of units throughout the paper in order to make the corresponding 
nonrelativistic limits as explicit as possible. The fundamental constants are speed of light in vacuum 
c = 2.99792458 x 10^° cms'^, Planck's constant h = 6.6260755 x 10"^^ J s, mass m = = 
9.1093897 X 10"^^ gm and electric charge e = \e\ = 1.60217733 x 10"^^ C of the electron, Bohr 
radius ao = fi^/me^ = 0.529177249 x 10~^ cm, Sommerfeld's fine structure constant a = /he = 
7.29735308 x 10-^ and Compton wave length X/27i = h/mc = 2.42631058 x 10"^° cm. 

We dedicate this paper to the memory of Professors A. F. Nikiforov and V. B. Uvarov in a 
hope that those masters would appreciate our effort to make their method more available for the 
beginners. 
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2.1. Evaluation of an Integral. Let us compute the following integral 

J:L= / e-^x-^^Ll{x)Li{x) dx, (2.1) 







where n > m and a — (3 = 0, ±1, ±2, ... . Similar integrals were evaluated in [T7], [28] and [5T], see 
also references therein, but an important relation with the Hahn polynomials seems to be missing. 

It is convenient to assume at the beginning that parameter s takes some continuous values such 
that a + s > —1 for convergence of the integral. Using the Rodrigues formula for the Laguerre 
polynomials [55], [57], [71] 

(x) = ^e^x"" (x"+"e-^) , (2.2) 
see the proof in Section 5 of the present paper, and integrating by parts 



J:L = -J (a;"+"e"^y"^ (x^L^ (x)) rfx 







1 [{x---e-i-'' [x^Ll (X))) i; - 1 / (x-"e-^y"'^^ (xU^ (x))' dx 



1 



But, in view of ([83), 



^ ' (x"+"e-") (x^L^(x))^"^ dx. 







i^sjp /^^^W_ ^(^ + m + l) .^ (-m), . fc+,^(n) 
^"^^-^"^^^ - m!r(/3 + l) 4. A:! (/3 + 1), ^ ^^"^^ 



r(/3 + m + l)r(s + l) ^ (-m)Js + l) 



X 



fc+s— n 



m!r(/3 + l)r(s-n + l)^A;! (/3 + 1), (, - n + 1), 
and with the help of Euler's integral representation for the gamma function [3], [57] 

/■oo 

/ x"+^+'e-^' dx = r (a + A; + s + 1) = (a + s + 1);^ F (a + s + 1) , 
see also (18. lip below, one gets 

.a/3 ^ r(a + g + i)r(/3 + m + i)r(5 + i) 

^ ^ n! m! r(/5 + l)r(s-n + l) 

(— m, s + 1, a + s + 1 \ 
(2.4) 

See [To] or equation (18.11) below for the definition of the generalized hypergeometric series 3F2 (1) . 
Thomae's transformation (18.91) . see also [10] or [12], results in 



00 



JnL = I e-'^x^^' LI (x) Li (x) dx (2.5) 
r (a + s + 1) r (/3 + m + 1) r (s + 1) 



ml (n — m)! F (/3 + 1) F (s — n + m + 1) 
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-m, s + 1, f3 — a — s 
X 3F2 I 1 , n > m, 

p + 1, n — m + I 

where parameter s may take some integer values. This estabhshes a connection with the Hahn 
polynomials given by equation fl8.7p below; one can also rewrite this integral in terms of the dual 
Hahn polynomials |55j . 

2.2. Special Cases. Letting s = and a = (3 in (12. 5p results in the orthogonality relation for the 
Laguerre polynomials. Two special cases 

I, = J^^^ = r e-^x"+^ (L^ {x)f dx={a + 2n + 1) ^ + + ^) (2.6) 



and 



-^2 — Jn, n-1, 2 ~ y *^ ^X°'L"_i (x) ^ (x) rfx — —2 ^ ^^^'^ 



(2.7) 



are convenient for normalization of the wave functions of the discrete spectra in the nonrelativistic 
and relativistic Coulomb problems [I?], [FFj . 

Two other special cases of a particular interest in this paper are 

POO 

Jk = JZk = / e--^-^' {L- {^)f dx (2.8) 



_ T{a + k + l)T{a + n + l) ^ / -k, k + 1, -n 
n! r(a + l) ' ' 1, a + 1 

and 

/"OO 

J-k-1 = JZ = / e-^x-^-' (L^ {x)f dx (2.9) 
Jo 

_T{a-k)T{a + n + l) ^ f -k, k + 1, -n 
~ n! r (a + 1) ' ' 1^ 1, a + 1 

The Chebyshev polynomials of a discrete variable tk (x) are special case of the Hahn polynomials 
4 (x, N) = hf °^ (x, A^) [75], [76] and [77|. Thus from (I231)-(E9]) and (IH2D one finally gets 



Jk = JZk= / e-^x^^' {Ll{x)y dx (2.10) 
Tia + n+l] 



tk (n, -a) 



OO 



and 

J_._i = JZ-k^i= I e-^x--'-' {L:{x)r dx (2.11) 

^0 

r (a — /c) r (a + n + 1) 

= TFT — I L I IN ^fc 

n! 1 (a + K + 1) 

for < A; < a. One can see that the positivity of these integrals is related to a nonstandard 
orthogonality relation for the corresponding Chebyshev polynomials of a discrete variable tk {x, N) 
when the parameter takes negative integer values N = —a. Indeed, according to the method of [55] 



THE HAHN POLYNOMIALS IN THE COULOMB PROBLEMS 



5 



and [57], these polynomials are orthogonal with the discrete uniform distribution on the interval 
[—a, —1] which contains all their zeros and, therefore, they are positive for all nonnegative values 
of their argument. The explicit representation (18.71) gives also a positive sum for all positive x and 
negative N. 



2.3. Connection Relation and Linearization Formula. The connection relation [3], [1] 

K = x: ~ (2.12) 

^-^ (n — m)\ 

is an easy consequence of the integral (12. 5p . we leave the details to the reader. 
A linearization formula 

n+m 

LI (x) [x] = J2 Cnmp («' 7) L; [x) (2.13) 

p=0 

gives a product of two Laguerre polynomials as a linear combination of other polynomials of the 
same kind; see |3], [1], [7], [61], [62] and references therein for a review on the linearization of prod- 
ucts of classical orthogonal polynomials. In view of the orthogonality relation, the corresponding 
linearization coefficients are given by 

^^^ + P + ^^ c„^^ (a, /3, ^) = r LI (x) Li (x) (x) x^e"^ dx, (2.14) 
V- Jo 

where one can assume that n > m and use the expansion 



Ll ix) 



p 



, . ^TttP^^' (2.15) 
P!r(7 + 1) t^/^!(7 + l), ^ ' 



by (18.21) . Then the integral of the product of three Laguerre polynomials is 



Cnmp (Q^) /3, 7) p (/y _|_ ]^) ^ ^ (/y _|_ ]^) ^ ^ "^^^ ^"^-^ ^"^-^ 

P 



= ^ V (~^^fe 7°^ r2 16) 

r(7 + l) ^A:!(7 + l), ^ • ^ 

and the remaining integral can be evaluated in terms of the Hahn or dual Hahn polynomials with 
the aid of (12. 5p . Positivity of the linearization coefficients is related to the orthogonality property 
of these polynomials. 

On the other hand, in the most important special case, 

n+m 

L:{x)LI{x)= Y1 <^nmpL;ix) (2.17) 

p=n—m>0 

with a = /? = 7, the linearization coeffitients Cnmp («) = Cnmp can be found as a single sum 

_( .xp (-p)n-m(« + l)m (-p)fc (-"^)fc T (2A; + 1) 

""^^ ^ ^ p\m\ 2^ k\{a + l),r{n-m-p + 2k+iy ^ 
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The summation is to be taken over all integer values of k such that 0<p — n + m<2k< 
min {2p, 2m) . 

Indeed, substituting ([23]) into I^AG^i 

/ _^ r(Q; + m + 1) >A (-P) 
r(a + l) Cnmp = (-1) — pT > 

ml In — mil ' ' I I « — n -U 

min(m,s) 

/c! (a + l)Jn-m + l) 



m! (n — m)! F (s — n + m + 1) 

fci -n V r-n - m -^ 1 V I • ^ 



with a = /3 = 7, then replacing the index s = n — m + /,0</<p — n + m and changing the order 
of summation one gets 

.n-m r (a + m + 1) 



r (a + 1) c„^p = (-1) — — (-p), 

m\ [n — m ! 



^ {-p + n - m)i{n - m + k + 1); ^qn 
^ /! (n-m- A; + 1); ^ ' ' 



with the help of 



{n — m + l)i 



(ri-m + / + l)^ = + 



— m — A; + 1); 



(n — m + l)j 

The 2-^1 (1) function in (12.201) is evaluated by a limining case of the Gauss summation formula (18.101) 

as 

—p + n — m, n — m + k + 1 
:-n + m),,f,| + ! ;1| (2.21) 



= (-ir'"-' (-p) 

This results in (12.181) and our proof is complete 



n-m-p^_^^, (n-m)! r(2A; + l) 

p\ T {n - m - p + 2k + 1)' 



Equations (I2.14p and (12.181) imply the following positivity property 



oo 



(_l)«.+n+p / ^^^-^ > ^2.22) 

when a > — 1 and m,n,p = 0, 1,2, ... ; see problem 33 on p. 400 of [3j. 

The single sum in f l2.18p can be rewritten as a 4F3. There are distinct representations for even 
and odd values ofe = (p — n + m)/2. When e = the summation formula (18.101) gives 

(n,\ - r(a + n + l) 

ml T [a + n — m + 1) 
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and the case n = m corresponds to correct value of the squared norm of the Laguerre polynomials. 
If e = {p — n + m) /2 is a positive even number the result is 

-P)n-m i-P)e i'^) e F (26 + 1) T (a + m + 1) 

p\ ml r (e) r (a + e + 1) 

-p + e, -m + e, e + 1/2, e + 1 \ 
X4i^3l 1/2, a + e + l, e 1 (2-24) 

(-P)n-m i-P)e+l/2 (-"^).+i/2 T (2e + 2) T (tt + m + 1) 

^ ' pi ml r (e + 1/2) r (a + e + 3/2) 

/ -p + e + 1/2, -m + e + 1/2, e + 1, e + 3/2 \ 
X 4^3 3/2, a + e + 3/2, e + 1/2 ' (2-25) 



and 



Cnmp (c^) 



if e = (j9 — 72 + m) /2 is an odd. 

Although we have not been able to find the linearization formula for Laguerre polynomials in 
explicit form in the literature, a closed expression for the linearization coefficients should follow as 
a limiting case of equations (1.7)-(1.8) for the linearization coefficients of the Jacobi polynomials in 
Rahman's paper [61]; see also [62] for g-extension of his result. Dick Askey has told us that he had 
computed the corresponding integral of the product of three Laguerre polynomials with the help 
of the generating function for these polynomials as a different single sum. In his opinion, this may 
have been found in the 1930's by Erdelyi or someone else. 

2.4. An Extension. The integral (12.11) has a somewhat useful extension 

JuL (^) = / e--a;°+^L° (x) (x) dx, (2.26) 



where n > m and a — [3 = 0, ±1,±2, ... and J^msi^) ~ Jnms- This integral is evaluated in the 
following fashion. Upon integrating by parts 

J:L{z) = e-^(-l)^^^^."+%"^,n.)(.U^(.))(^"^) (2.27) 

k=l 

^ r (x"+"e--y"-') (x^L^(x))^'^ dx 

J z 

and putting / = n we use (12.31) and the integral representation for incomplete gamma function [36] . 

T{a,z)= e-^'x"-^ dx, Re a > 0. (2.28) 

J z 

As a result 

(^) = / e--a;"+U;^ (x) (x) dx (2.29) 

Z [Zj \Z 1^-^ \Zj , 



k=l 
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n] m\ ^—^ k! (p + 1)^ i [s — n + k + 1) 



k= 

with n > m and Re {a + s) > —1. Here one can use the standard relations [36], [57] : 

T{a,z) = e^'G{l- a,l- a,z) (2.30) 

= T{a) -- z"F{a,l + a,-z) 
a 

= r (a) - ^" V , ,~ ^ ^, , , Rea>0 

with the confluent hypergeometric functions F(a,7, and G{a,'y,z) , respectively. The last ex- 
pression gives the asymptotic of the incomplete gamma function as z — and |arg(2;)| < vr. The 
asymptotic for large values of z is 

T{a,z) = e-' z'^G {1,1 + a, z) (2.31) 

^ .-^.-f£(-i)'(izi!).^o(i 

\fc=0 ^ 

as z ^ oo and |arg(z)| < vr; see [3B] and ^7\. This allows to find asymptotic expansions of the 
integral fl2.29p as 2 ^ and z ^ 00. Integration by parts in fl2.28p results in the functional relation 

r(a + l,2) = a^(a,z) + z"~^e-^ (2.32) 

which implies 

n ^ 

T{n+l,z)=n\e-^Y.'k\ ^^"^^^ 

k=0 

for positive integers n. 

3. NONRELATIVISTIC COLOUMB PROBLEM 

3.1. Wave Functions and Energy Levels. The nonrelativistic Coulomb wave functions obtained 
by the method of separation of the variables in spherical coordinates, see Section 7.1, are 

i^=i^nlmi^) = Rndr) Yi^{9,^), (3.1) 

where Yim {0, (p) are the spherical harmonics, the radial functions Rni (r) are given in terms the 
Laguerre polynomials [17], [51], [57], [61] 



with 

2Z / r\ 



ao = — - (3.3) 



and the normalization is 



00 



i?^,(r)r^rfr = l. (3.4) 







Here n = 1, 2, 3, ... is the principal quantum number of the hydrogen-like atom in the nonrelativistic 
Schrodinger theory; / = 0, 1, ... , n — 1 and m = —I, — / + 1, ... ,1 — 1,1 are the quantum numbers of 
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the angular momentum and its projection on the z-axis, respectively. The corresponding discrete 
energy levels in the cgs units are given by Bohr's formula 

where n = 1,2, 3, ... is the principal quantum number; they do not depend on the quantum number 
of the orbital angular momenta / due to a "hidden" SO (4) —symmetry of the Hamiltonian of the 
nonrelativistic hydrogen atom; see, for example, [M], [51] and references therein and the original 
paper by Fock [H] and Bargmann [11] . 



3.2. Matrix Elements. In this section we evaluate the mean values 

|2 



r'P) = ^ = , dv = r^drduj (3.6) 

J Jo 

in terms of the Chebyshev polynomials of a discrete variable tk {x, N) = h^^' {x, N) [75] , [76] and 
[77]. Here we have used the orthogonality relation for the spherical harmonics ^7], |80j . 



{6, ip) Yi,^, {e, if) duj = 6w6mm' (3.7) 
with du = sin 6 dOdip and O<0<7r,O<(y9< 27r. The end results are 

(r'=-^> = ^(^)'"t.(n-/-l,-2/-l), (3.8) 



when k = 0, 1, 2, ... and 



when A; = 0, 1, 21. 



1 / 2Z ^ ''^^ 



1 2n \naQ 



tk{n-l- 1,-21-1), (3.9) 



Although a connection of the mean values (13. 6p with a family of the hypergeometric polynomials 
was established by Pasternack [58], the relation with the Chebyshev polynomials of a discrete 
variable was missing. This is a curious but fruitful case of a "mistaken identity" in the theory 
of classical orthogonal polynomials. The so-called Hahn polynomials of a discrete variable were 
originally introduced by Chebyshev [77j, they have a discrete measure on the finite equidistant set 
of points. Bateman, in a series of papers [12], [13], [H], [15], and Hardy [H] were the first who 
studied a continuous measure on the entire real line for the simplest special case of these polynomials 
of Chebyshev. Pasternack gave an extension of the results of Bateman to a one parameter family 
of the continuous orthogonal polynomials [59]. After investigation of these Bateman-Pasternack 
polynomials in the fifties by several authors; see [79], [86], [23], [25] and [26]; Askey and Wilson [6] 
introduced what nowadays known as the symmetric continuous Hahn polynomials, they have two 
free parameters — but one parameter had been yet missing! Finally, Suslov [HZ], Atakishiyev and 
Suslov [S] and Askey [5] have introduced the continuous Hahn polynomials in their full generality 
in the mid of eighties. More details on the discovery the continuous Hahn polynomials and their 
properties are given in [39] among other things. 
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Indeed, in view of the normalization condition of the Coulomb wave functions (13 ■4p one gets 

oo 

(rP) = [ Rli{r)rP+^ dr (3.10) 



and the last integral can be evaluated with the help of fl2.10p or (12. Ill) giving rise to fl3.8p and fl3.9p . 
respectively. 



A convenient "inversion" relation for the Coulomb matrix elements 

2^x2fc+i (2/ -A;)! 



r^+^/ ynaoj {21 + k + 1)1 ^ ^ ^ ' 

with < A; < 2/, follows directly from (13. 8p and (13. 9p . This relation is contained in an implicit form 
in [51], it was given explicitly in |58j . 

3.3. Special Cases. The explicit expression (13. Sp for the matrix elements (r^) and the familiar 
three term recurrence relation for the Hahn polynomials h^j^' (x, N) [55], |57j . 

xhi"^ ^) (x, N) = a,h^^;f^ (x, N) + P,hi"^ (x, N) + 7,^!:/) (x, N) (3.12) 



with 



{n+l){a + /3 + n + l] 

Oik =- 



{a + (3 + 2n + 1) {a + (3 + 2n + 2)' 
a-(3 + 2N -2 {0" -a^){a + ^ + 2N) 



4 A{a + f3 + 2n) {a + (3 + 2n + 2)' 

{a + n){(3 + n){a + p + N + n){N -n) 
{a + P + 2n){a + P + 2n + l) ' 

imply the following three term recurrence relation for the matrix elements (13. 6p : 



Ik 



+ 

\2 



k{i2l + lY-e) fnaoy^^,^,^ 



k + 1 \2Z 
with the "initial conditions" 



-) = A' (1) = 1 (3-14) 

which is convenient for evaluation of the mean values (r^) for k > 1 [58]. The inversion relation 
(13. lip can be used then for all possible negative values of k. One can easily find the following matrix 
elements 

(r) = ^(3n2 -/(/ + !)), (3.15) 
(r^> = 2(^)'(5n^ + l-3/(/ + l)), (3.16) 
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^ ^ " (3.17) 



r / aov? 



r2 / agn3 (2/ + 1) ' 



2^2 ^ ^ 

(3.18) 



rV a^n^ (/ + !)(/ + 1/2)/^ 
1 \ (3^2 - / (/ + 1)) 



r.4 



(3.19) 

2aX + 3/2) (/ + 1) (/ + 1/2) / (/ ^T72) ' ^^'^^^ 
which are important in many calculations in quantum mechanics and quantum electrodynamics [1], 
[T7] . [16], [17], [82]; see [17] for more examples. 

Equations (13.51) and (13.171) show that the total energy of the electron in the hydrogen-like atom 
is equal to half the average potential energy: 

/1\ Z2p2 

{U) = -Ze^ (-) = = 2E. (3.21) 

\ r / ao^ 

This is the statement of so-called virial theorem in nonrelativistic quantum mechanics; see, for 
example, [T7], p. 165 and [51] . 

The average distance between the electron and the nucleus r = (r) is given by (I3.15p . The mean 
square deviation of the nucleus-electron separation is 

= ^ - = (S)' + 2) -1^1 + if) . (3.22) 

The quantum mechanical analogue to Bohr orbits of large eccentricity corresponds to large values 
of this number (small I) . 

3.4. Screening. Let us evaluate the effective electrostatic potential V (r) created by motion of 
the electron in a hydrogen-like atom with the nucleus of charge Ze. This result is well known in 
the nonrelativistic Schrodinger theory — see, for example, [l5] and [5^ — but we emphasize the 
connection with the Hahn polynomials in order to obtain similar results in the relativistic Dirac 
theory in the next section. For the electron in the stationary state with the wave function (13. ip and 
the quantum numbers n, I and m this potential is 

V{r) = — -e! [r'f dr'du', (3.23) 



where ep (r) = e lipnim is the average charge distribution of the electron in the atom. In order 
to evaluate the integral one can use the generating relation (18.150 in the form 



°° S / A ^ 



^ I 2s + 1 ^ 

m'=—s 



Y:^,{e',v')Y,^,{e,^)], (3.24) 

=0 '> ' ■ '^'"^ ' " 

where r< = min (r,r') and r> = max (r,r'); see [57], [63], [80] for the proof of this expansion 
formula. In view of (13. ip one gets 

/ ('-''' = 5: 5^ / 5t flS, (/) (r' f dr' (3.25) 



R3 
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X J2 Ysn^'iO,^) I Y:^,{9',^')Y:^{9',^')Y,mie',^') dco'. 

m'=—s 

The (fi' integration reduces the sum over m' to a single term with m' = and we arrive at 

oo °° 



x^o {0, 0) j F/o YL (0', Yim {e\ dcu'. (3.26) 



S2 



The integral of the product of three spherical harmonics can be evaluated in terms of the Clebsch- 
Gordan coefficients C/™ , of the quantum theory of angular momentum [35], [55], [55], [63], [80] , 



■'Iimil2m2 

[8T] . [8^ with the help of the product formula (18.161) and the orthogonality property (13.71) . The 
result is 



Yso iO', YL {0', Yim {e\ du' = CiZo Qo.o- (3.27) 

In view of the symmetry property of the Clebsch-Gordan coefficients [55], [63], [80] 

and the selection rule |/i — /2I < ^ < ^1 + ^2 of the addition of two angular momenta in quantum 
mechanics, the integral (I3.27P is not zero only for s = 0, 2, ... , 21. As a result 

I 



00 

X / ^,Rli{r'){r'fdr'. (3.29) 


It is known that the Clebsch-Gordan coefficients are simply the Hahn polynomials up to a nor- 
malization factor; see [55] and references therein for more details on the relation between the 
Clebsch-Gordan coefficients and the Hahn polynomials, which was overlooked on the early stage 
of developing of the quantum theory of angular momentum [27], [2S]) [S3], [SD], [SI], [S3], [SI] and 
had been established much later independently by Koornwinder [SO] and Smorodinskii and Suslov 
[65] . It is worth noting that before that Wilson [85] found that the next "building blocks" of the 
quantum theory of angular momentum, the so-called 6j-symbols, are orthogonal polynomials of 
a discrete variable, see also [66], [68] and [69]; Askey and Wilson ^ studied this new orthogonal 
polynomials and their g-extensions in details; see also [2], [3], [12], [3S], [55] and references therein 
for the current status of the theory of Askey- Wilson polynomials and their special and/or limiting 
cases. Thus the Hahn polynomials appear in the expression fl3.29p for the effective electrostatic 
potential V (r) in the nonrelativistic hydrogen-like atom. 

The special Clebsch-Gordan coefficients C/q 2^ (13.291) are 



.10 ,,sil + s)\i2s)\ I (21 + 1) {21 - 2s)\ 



(/-.)!(.!)^V (2/ + 2. + 1)! ^^-^^^ 
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and 

y2s, {e, ^) = P2S (cos e) , (3.31) 

where P„ (x) are the Legendre polynomials. 

The integral over the radial functions in (13.291) can be rewritten in the form 



oo 

„2, 



RlAr') {r'f dr' (3.32) 



Jo Jr 



{rT^'RlAr')dr'- I {r'^^' Rl,{r') dr' 



J.2S+1 







/■oo 

+r2^ / {r'f-''Rl{r')dr' 

J r 



^ - (r) + r^^ J_2.-i (r) 



where the first integral = (^^^) is given by (13. 8p in terms of the Chebyshev polynomials of a 
discrete variable. The other two integral are of the form 

POO 

Jk{r)= {r'f^'Rl,{r')dr'. (3.33) 



They are special cases of our integral (I2.26P and can be evaluated by (12.291) in terms of the incomplete 
gamma function; they have simple asymptotics at infinity. The reader can work out the details. 

For the electron in the ground state n = 1 and / = m = all the integral are easily evaluated 
and the result is 

V (r) = +('- + ^\ e--^/-. (3.34) 



r \r aQ 

For small distances r ^ the effective potential V {r) —>■ eZ/r as expected and as r ^ oo the 
potential V (r) — > e{Z — 1) /r which is the potential of the nucleus of charge Ze screened by the 
electron. 

4. Relativistic Coulomb Problem 
4.1. Dirac Equation. The relativistic wave equation of Dirac [31], [32], [33], [38], [39] 



d 

ih—ip = Hi: (4.1) 

for the electron in an external central field with the potential energy U (r) has the Hamiltonian of 
the form 

H = ccxp + mc^p + t/ (r) , (4.2) 
where ap = aipi + 0:2^2 + «3P3 with the momentum operator p = — i/iV and 

-ill)' ^-il-.): ^-il)- («) 
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We use the standard representation of the Pauh matrices 



and 





1 






< -i 






w 


(72 = 1 


y i 



0-3 



(4.4) 












1 
1 



The relativistic electron has a four component wave function 

/ ^i(r,t) \ 



ip = ip{r,t) 



u(r,t) 
v(r,t) 



(4.5) 



■^2 (r,t) 
^3 (r,t) 
V V^4(r,t) / 

and the Dirac equation (14. ip is a matrix equation that is equivalent to a system of four first order 
partial differential equations. The inner product for two Dirac (bispinor) wave functions 

V'2 
V^4 / 



Vl 



"2 
V2 



02 
<^3 

04 / 



is defined as a scalar quantity 
with the squared norm 



/ V'V dv = (ui"''u2 + vi"''v2) dv 

^101 + ^202 +^3 03 +^4 04 



(4.6) 



R3 



dv 



/ ip'^ip dv = I (ui''"ui + vi''"vi) dv 

JR3 JR3 



(4.7) 



R3 



+ l^2l' + IV^3l' + IV'4l') dv 



and the wave functions are usually normalized so that 
The substitution 

gives the stationary Dirac equation 

(r) = E^Ij (r) , 

where E is the total energy of the electron. 



1/2 



1. 



(4.8) 
(4.9) 



According to Steven Weinberg ([82j, vol. 1, p. 565), physicists learn in kindergarten how to solve 
problems related to the wave equation of Dirac in the presence of external fields. In Section 6 of 
this paper, for the benefits of the reader who is not an expert in theoretical physics, we outline 
a procedure of separation of the variables and solve the corresponding first order system of radial 
equations for the Dirac equation in the Coulomb field U [r) = —Ze^/r. The end results are presented 
in the next section; see also [1], [I7j, p^, ^39j, [47], [5l], [57], [M] and references therein for more 
information. 
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4.2. Relativistic Coulomb Wave Functions and Discrete Energy Levels. The exact solu- 
tions of the stationary Dirac equation 

= Eip (4.10) 

for the Coulomb potential can be obtained in the spherical coordinates a result of a rather lengthy 
calculation. Remarkably the energy levels fl4.16p were discovered in 1916 by Sommerfeld from the 
"old" Bohr quantum theory and the corresponding (bispinor) Dirac wave functions were originally 
found by Darwin [30] and Gordon [43] at early age of discovery of the "new" wave mechanics; see 
also [T8] for a modern discussion of "Sommerfeld's puzzle". These classical results are nowadays 
included in all textbooks on relativistic quantum mechanics, quantum field theory and advanced 
texts on mathematical physics; see, for example, pQ, [TB], [I?], [22], [IZ], [SI] and [SZ]- More details 
are given in the last but one section of this paper. The end result is 

3^t(n) F{r) 



X 



iyj^ (n) 



G{r) 



(4.11) 



where the spinor spherical harmonics yf^ (n) = y 



ordinary spherical functions Yim (n) , n = n (6^, yj) 
with the spin 1/2 as follows [T], [IS, [63], [80 

/ 



(i±i/2) 

' jr'm 



yfm (n) 




(j + 1/2) T{m- 1/2) 
2j + (1 ± 1) 

(j + 1/2) ± (m + 1/2) 
2j + (1 ± 1) 



n) are given explicitly in terms of the 
r/r and the special Clebsch-Gordan coefficients 



^■±1/2, m-1/2 (n) 



^■±1/2, m+1/2 (n) 



with the total angular momentum j = 1/2,3/2,5/2, ... and m = + 1, 

Section 6.1 for discussion of properties of the spinor spherical harmonics in detail. 



(4.12) 



,j - l,j; see 



The radial functions F (r) and G (r) can be presented as [57] 




(ek — v) n\ 



Here (^) are the Laguerre polynomials and we use the following notations k 



fi 
9i 




V 



a/k^ — /i^, /i = (xZ = Ze'^/hc, a = Vl — e^, e 

^ = 2apr 



E/mc , (3 = mc/h = 27r/A and 



2Vl - — r. 
n 



The elements of 2 x 2-transition matrix in (14.131) are given by 

a/i „ a(K — u] 



fi 



f2 



9i 



92 = 



(4.13) 

± (j + 1/2) , 

(4.14) 
(4.15) 



en — u SK — u 

This particular form of the relativistic radial functions (14. 13p is due to Nikiforov and Uvarov |57j ; 
it is very convenient for taking the nonrelativistic limit c —>■ oo, see more details in Section 6.4; we 
shall use this form throughout the paper instead of the traditional one given in [T], [TB], QZj, ^28] 
and elsewhere; cf. equation (16.821) below. 
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The relativistic discrete energy levels e = En = En/ Eq with the rest mass energy Eq = mc^ are 
given by the famous Sommerfeld fine structure formula 



En = I iJ, = aZ = — — , = v^K^ - (4.16) 

^1 + /iV {n + vf 

Here n = = 0, 1, 2, ... is the radial quantum number and k = ± (j + 1/2) = ±1, ±2, ±3, ... . 
In the nonrelativistic limit c ^ oo one can expand the exact Sommerfeld-Dirac formula (I4.16p in 
ascending powers of /i^ = (aZ)^ , the first terms in this expansion are 

E ^ f n 3^ 



mc 



2 2n2 2n^ 



Here = ra^ + j + 1/2 is the principal quantum number of the nonrelativistic hydrogen-like atom. 
The first term in this expansion is simply the rest mass energy = of the electron, the second 
term coincides with the energy eigenvalue in the nonrelativistic Schrodinger theory (13.51) and the 
third term gives the so-called fine structure of the energy levels - the correction obtained for the 
energy in the Pauli approximation which includes interaction of the spin of the electron with its 
orbital angular momentum; see |T7] and [M] for further discussion of the hydrogen-like energy levels 
including comparison with the experimental data. One can show that in the same limit /i — > the 
relativistic Coulomb wave functions (14. lip tend to the nonrelativistic wave functions of the Pauli 
theory; see, for example, [SZ] for more details; we shall elaborate more on this limit in Section 6.4. 

We give below the explicit form of the radial wave functions (14.131) for the 15'i/2-state, when 
n = = 0, / = 0, j = 1/2 and n = —1 : 

/ -1 \ 



E{r) \ /2Zy/' I z/i + 1 



G(r) \ao J V 2r (2z/i + 1) 



1-z/i 



V V 1 + ^1 / 



Here ui = ^Jl — n'^ = Si and = 2^yl — e\l3r = 2Z {r/ao) . One can see also [T], [T7], [16], [30] , 
[29] . [43] . [17], [5l], [SI] and references therein for more information on the relativistic Coulomb 
problem. 

4.3. Matrix Elements. In this section we evaluate the mean values 



(rO = ^ = (4.19) 



3 







of all possible powers of r with respect to the relativistic Coulomb functions given by (I4.1ip - P:.13p in 
terms of the Hahn polynomials (18. 7p . First we use the orthogonality relation of the spinor spherical 
harmonics; 

{y^i (n))^ y^Z (n) duj = 5,,,5a'5mm', (4.20) 
see [1] , [16] , [80] or (16.60 ; and the explicit form of the wave functions (14. lip in order to simplify 
j ^'^ rPtl^ rfy = y (3;t (n) y (n)) {F^ (r) + (r)) r^drdiu 

R3 R3 



THE HAHN POLYNOMIALS IN THE COULOMB PROBLEMS 17 

//'CO 
{y^{n)y{n)) duj j (r) + (r)) dr 



52 



r-P+^ {F^ {r) + G'' {r)) dr. (4.21) 

More details on construction of the spinor spherical harmonics 3^j^ (n) = D^j^^^^"* (n) and study 
of their properties are given in Section 6.4. The radial wave functions in f l4.13p are normalized as 
follows 

/•oo 

/ r\F\r) + G\r)) dr = l. (4.22) 
Jo 

Thus one needs to evaluate the integral 



00 



rP) = / rP+^ [F^ (r) + (r)) dr (4.23) 



only, and our final result with the notations from the previous section can be presented in the 
following closed form 

Aiip'^{2a[5Y {r'P) (4.24) 

,r(2z/ + p + 3) / 1-n, p + 2, 

= an [en + z/ — — r — -xFo 

^ ^ r(2z/ + 2) ' ' \^ 2Z/ + 2, 1 

, o^ / 2 2 2Nr(2z/ + p + 2) /l-n, p + 3, 

-2 (n + 2) u - v'^) —— — 3F2 

> r(2z/ + 2) ^ 2Z/ + 2, 2 

,r(2z/ + j9+l) /-n, p + 2, 

+ aK £K — Z/ ; -tPo 

^ ' T{2u) ' ' 1^ 2z/, 1 

where the terminating generalized hypergeometric series 3F2 (1) are related to the Hahn and Cheby- 
shev polynomials of a discrete variable due to fl8.7p . 

Substituting Km into KTS^ one gets 

/oo 
^_C^2.+p (^) + f^Ll^-i (4 25) 

+ (^7ie«MO+^?2^r^(0)') 

POO 

^0 

/■oo 

I'OO 



where 



C = 8^z/2 20/?^ ^ / \ , 4.26 

a (EK — V) n\ 
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and, in view of (14.151) . 

fi+9l = —^ 72^, hf2 + 9ig2 = fi + gl = 2K{K-v). (4.27) 

[SK — uj en — u 

The two types of the integrals occuring in the calculation are given by our "master" formula (12.51) 
as follows 

I e-^e''^-^' (4.28) 


T{2u + p + 3)T{2u + n + l) ( 1 - p + 2, -p-l 



(n-l)\V{2u + 2) ' M 2Z/ + 2, 1 



Jo 



(4.29) 



(p + 2) r (2z/ + p + 2) r (2z/ + n + 1) / 1-n, p + 3, -p 

3-^2 



(n- l)!r(2z/ + 2) " M 2Z/ + 2, 2 



and 



/■oo 

/ e-^e-'-'HLn-HOT (4.30) 
^0 



r (2z/ + p + 1) r (2z/ + n) / -n, p + 2, -p-l 



n\T {2v) ' ' \^ 2z/, 1 

Substituting these integrals into (I4.25P and using the identity 

c?n{2v^n)=e^K^ -v^, (4.31) 
see Section 6.3 for the proof, one finally arrives at (I4.24p . 

In view of (18.71) . the relations with the Hahn and Chebyshev polynomials of a discrete variable 



are 



and 



4/ii/2 {2a^Y {r^) = i^^^ + ^) (^ - 1, -1 - 2v) (4.32) 

-2^-±l^ (e\^ - /,{!, 1) _ 1 1 _ 2u) 

p+r^ / p V ' ^ 

+aK {sK — u) /ip^ili ^ (ri, 1 — 2v) , p > 
/ U=aK(.K + .)-i-PpZL,4%»>(„- 1.-1-2.) (4.33) 



{2apf^^ \rP+3/ ^ ^r(2z/ + p + 3) 



The averages of for the relativistic hydrogen atom were evaluated by Davis ^28j in a form which 
is slightly different from our equation (I4.24p . It does not appear to have been noticed that the 
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corresponding 3F2-functions can be expressed in terms of Hahn polynomials. The ease of handling 
of these matrix elements for the discrete levels is greatly increased if use is made of the known 
properties of these polynomials [36j, [IH], [55], [57] and [75], [76], [77]. For instance, the difference- 
differentiation formula 

A/i(^' P){x,N) = {a + (3 + m + 1) ^^'^ (x, iV - 1) , (4.34) 

where A/ (x) = / (x + 1) — / (x) , in the form 

/.(%°) (n, -2Z/) - /.(°;°) (n - 1, -2u) = {p + 2) h^^^ {n - 1, -1 - 2v) (4.35) 

allows to rewrite formulas fl4.32l) - fl4.34l) in terms of the Chebyshev polynomials of a discrete variable 
tm {x) = hm (x, A^) only. Use of the three term recurrence relation f l3.12p simplifies evaluation 
of the special cases of these averages. Equation (18.81) gives asymptotic formulas for the matrix 
elements as oo. 

4.4. Nonrelativistic Limit. In the limit c — oo the relativistic Coulomb matrix elements given 
by (I4.19P and (14.241) tend to the nonrelativistic ones (13. 6p . This can be easily shown with the help 
of the asymptotic formulas 



eK±u=iK±\K\)- '7,/" ,7" f^' + O (/) (4.36) 
as jU = Ze^ /he 0; see Section 6.4 for more details. 

4.5. Special Cases. Some important special cases of the relativistic matrix elements are 

/ 2\ ^ 5n (n + 2u) +Au^ + l-eK {2eK + 3) 

' 2 [a^f ^ ■ ^ 

(r) = ^ (3en {n + 2v) + k {2eK - 1)) , (4.38) 

(1) = 1, (4.39) 

= A(l-.^)(.z. + /.v/r^), (4.40) 

1 \ _ 2c? (5" K C2eK - 1) 

rV /iz/(4z/2-l) ' ^ - ' 

lUaCa/^f ^^X-^tr:'^ . (4.42) 



I ' z/ (z/2 - 1) (4z/2 - 1) 
Derivation of these closed forms requires more work than in the nonrelativistic case because of a 
more complicated structure of the general expression (I4.24p . Here are some calculation details. 

The special case p = implies the normalization of the radial wave functions (I4.22p and (I4.39p . 
Indeed, 

f I —n 2 — 1 \ 

4/iz/2 (1) = -4/i {e^K^ - 1/2) + an (e^ + v) {2v + 2) 3F2 ' ' (4.43) 



2z/ + 2, 1 



+an {en — v) 2v 3F2 



-n, 2, -1 
2z/, 1 
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= — 4/i (e^K^ — z/^) + 2aK [en + u){v + n) + 2aK [sk — u) (u + n) 
= AeK^a{v + n) - Aii[e^K^ - i/) = 
in view of the quantization rule 

efi = a{u + n) , (4.44) 

which leads to the Sommerfeld-Dirac formula for the discrete energy levels (14.161) : see Section 6.3 
for more details on (I4.44p . 

The special case p = — 1 of (14.241) reads 

/1\ / x / \ , A 

an [EK + u) + an [ek — u) (4.45) 

, . . r(2z/ + 1) / 1 - n, 1 

^ r 2Z/ + 2 ' M 2Z/ + 2 ' 



aP \r 



where by the summation formula of Gauss (IS.lOp 



2ae/t^ - 2/ia^n, (4.47) 



Thus, in view of (I4.3ip 

2p^ /I 

ap \r, 

and the final use of the quantization rule (14.440 results in (I4.4ip . Our calculation of (r^^) shows 
that there is no simple form of the virial theorem in the Dirac theory of relativistic electron moving 
in the central field of the Coulomb potential. 

In the case p = — 2 we get 

which leads to fOB . 

In a similar fashion, for p = —3 : 

i</lU2,(.V-.^)^i^ (4.49) 
{2a(3f\r'/ ^ T {2u + 2) ^ ^ 

, r(2z/) ^ / 2, -1 

+^^^'^^''^IWTY)'^'[ 2U + 2, 1 

, r(2z/-2) ^ f -n, 2, -1 

+aK (EK — U) — ^ qFo 

^ ' r (2z/) ' 2v, 1 

aK {ek + v) {2v + 2n) ^ an {ek — u) {2v + 2n) ^ 2/i {e'^k^ — 



(2z/ + 2) (2z/ + 1) 2z/ 2z/(2z/- l)(2z/-2) {2v ^ \)2v {2v - \) 
which is simplified to (14. 42 p . 
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The special case p = 1 gives the average distance r = (r) between the electron and the nucleus 
in the relativistic hydrogen-like atom. One gets 



4/iz/2 {2ap) (r) = -6/i {e\^ - u^) {2u + 2) 3F2 
+aK + z/) (2z/ + 2)2 3^2 
+aK {sK - u) (2z/)2 3F2 



1-n, 4, -1 
2z/ + 2, 2 

1 -n, 3, -2 
2z/ + 2, 1 

-n, 3, —2 
2u, 1 



where 



1 _ ^ 4 _i 



-n 3 —2 



Thus 



(3n (n + 2z/) + 2//^) - anu^ - 3/i (e^ 2)(z/ + n), (4.50) 



which can be simplified to fl4.38p by a straightforward calculation with the aid of (14.311) and (I4.44p : 
we leave the details to the reader. 



For p = 2 : 



where 



and 



4,.= i2apf ^ -8, - £(|±| ( I (4.51, 

, r(2z/ + 5) ^ / 1-n, 4, -3 
+^^^'^ + '^lWTY)-^^^[ 2. + 2, 1 

, r (2z/ + 3) ^ / -n, 4, -3 

^ ^ r (2z/) ' M 2i., 1 



(2i/ + 2)2 3i^2 I 2 1 =5n2 + 10w + 4z/2 + l 



_^ 4 _3 

(21^)3 3F2 I ' ' 1 = 4 (n + z/) (5^2 + lOnz/ + 2//^ + 3i^ + l) 



2z/, 1 

This can be transformed to (14.371) with the help of (I4.44p . 
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4.6. Screening. Let us evaluate the effective electrostatic potential V (r) for the relativistic hydro- 
gen-like atom. For the electron in the stationary state with the wave functions (14. lip corresponding 
to the total angular momentum j, its projection m and the radial quantum number n = rir this 
potential is 

n^)-'i-^f (4.52) 

where 

ep (r) = (r) ^ (r) = eQjm (n) {F^ (r) + (r)) , (4.53) 

Q,„(n) = (3^±„(n))^ yf^in) 

is the charge distribution of the electron in the atom. We evaluate the integral with the help of the 
generating relation ( 13.24p . Indeed, 



OO 



R3 ^ 

J2 Ysm' (n) [ Y;^, (n') Q,^ (n') du'. (4.54) 



X 

m'=—s 



In view of ( KW) and flS^T]) . 



i-i/2 



p=o ^ P 

and with the help of the orthogonality property (13. 7p one gets 

r ■'^^^^ ; 4^ 

y O (n') Qjm (n') c/w' = 5„/o Yl Y 2p. 



52 



p=0 

Substitution to (I4.54p gives 



CO 

^2s 



^ (2J + 2S + 1) (2j-2.)!(j + .-l/2)!(2s)! 

' (2j + l)(2j + 2s)! (^■_s-l/2)!(s!)' 

{F'{r') + G'{r')){rfdr'. 

J > > 


with the help of fl638D . 

The integral over the radial functions in (14.560 can be rewritten in the form 



00 

' < 



' > 



(F' (r') + (/)) {rfdr' (4.57) 
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f\2s+2 



1 







y,2s + l 







J.2S+1 



(^')2^+2 (^2 ^ ^2 (^/)^ 

o 

+r2^ / (r')'"'' (F^ (r') + (r')) cir', 

where the first integral is given by fl4.24p . The next two can be evaluated with the help of (12.291) . 
For the electron in the lS'i/2-state with the radial functions (14.181) the result is 

V{r) = ^" + ^^~n'. 1 X r""^-^e-"^'"/°° (4.58) 

+ 



r r(2z/i + l) 

T {2ui,2Zr/ao) [ 2uie 2Ze 



r(2z/i + l) \ r ao 



where i^i = y 1 — jj,"^. For small distances r — the effective potential V (r) — >■ eZ/r and as r — oo 
the potential V (r) —* e{Z — 1) /r which is the potential of the nucleus of charge Ze screened by 
the electron. In the limit c ^ oo one gets the nonrelativistic formula (I3.34p . 

5. Special Functions and Quantum Mechanics 

In this section we give a short summary of Nikiforov and Uvarov's approach to special functions 
of mathematical physics and their applications in quantum mechanics j57j . 

5.1. Generalized Equation of Hypergeometric Type. The second order differential equation 
of the form 

«"+^«'+4tt'' = 0. (5-1) 

a [z) cr"^ [Z) 

where a (z) and a [z) are polynomials of degree at most 2 and r (2;) is a polynomial of degree at 
most 1 of a complex variable z, is called the generalized equation of hypergeometric type. By the 
substitution u = ip{z)y equation (15. ip can be reduced to the equation of hypergeometric type 

a{z)y" + T{z)y' + \y = Q, (5.2) 

where r [z) is a polynomial of degree at most 1, and A is a constant. The factor ip {z) here satisfies 

^ = ^, (5.3) 

V a (2;) 

where tt (2;) is a polynomial of degree at most 1 given by a quadratic formula 



T,{z) = ^—^±\li'^—^\ -a + ka (5.4) 
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and constant k is determined by the condition that the discriminant of the quadratic polynomial 
under the square root sign is zero. Then r (z) and A are determined by 

t{z) = t{z) + 27i{z) , \ = k + 7r'{z). (5.5) 

Two exceptions are 



(1) If a (z) has a double root, a (z) = [z — of' , the original equation can be carried out into a 
generalized equation of hypergeometric type with a (s) = s, by a substitution s = {z — a)'^ . 

(2) If a (yz) = 1 and (j {z) /2)^ — a [z] is a polynomial of degree 1, the substitution tt (z) = 
— r {z) /2 reduces the original equation to the form 

y" + {az + h)y = 0. (5.6) 

The linear transformation s = az + h takes this into a Lommel equation fl5.68p . 

Solutions of fl5.ip - fl5.2l) are known as special functions of hypergeometric type; they include 
classical orthogonal polynomials, hypergeometric and confluent hypergeometric functions, Hermite 
functions, Bessel functions and spherical harmonics. These functions are often called special func- 
tions of mathematical physics. 

5.2. Classical Orthogonal Polynomials. The Jacobi, Laguerre and Hermite polynomials are the 
simplest solutions of the equation of hypergeometric type. By differentiating (15.21) we verify that 
the function vi = y' {z) satisfy the equation of the same type 

a{z)vi +Ti{z)vi + ^i^Vi = Q, (5.7) 

where ri {z) = t {z) + o' (z) is a polynomial of degree at most 1, and = A + r' (z) is a constant. 

The converse is also true: any solution of (15. 7p is the derivative of a solution of (15.21) if A = 
yU]^ — r' 7^ 0. Let v\ (z) be a solution of (15. 7p and define the function 

y (z) = (a (2;) v'^^r {z) v^) ■ 

Then 

\y' = - (av'l + Tiv[ + t'vi) = \Vi 
or Vi = y' {z) and, therefore, y {z) satisfy (15.20 . 

By differentiating (15.20 n times we obtain an equation of hypergeometric type for the function 

a (z) Vn" + Tn {z) + /U^t^n = 0, (5.8) 

where 

Tn {z) =t{z)+ na' {z) , (5.9) 

/i„ = A + nr' + ^n(n-l)a". (5.10) 

This property lets us construct the simplest solutions of (15.21) corresponding to some values of A. 
Indeed, when /i„ = equation (15.80 has the trivial solution t>„ = constant. Since Vn [z^ 
the equation (15. 2p has a particular solution y = yn (z) which is a polynomial of degree n if 

X = Xn = -nr' --n{n-l)a" (n = 0, 1, 2, ... ) . (5.11) 
2 
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To find these polynomials explicitly let us rewrite equations (I5.2p and (IS.SI) in the self-adjoint 
forms 

(apy')' + Apy = 0, (5.12) 

(f^PnO' + ^J'nPnVn = 0. (5.13) 

Functions p (z) and p„ (z) satisfy the first order differential equations 

(ap)' = TP, (5.14) 

i^Pn)' = '^nPn- (5-15) 



So, 



whence 



and, consequently, 

p„ (^) = a" (^) p (z) . (5.16) 
Since crp„ = p^^;^ and = Vn+i one can rewrite (15.131) in the form 



{p^+iVn+l)' . 



Hence we obtain successively 



Pn^n 



py = PqVq = -— (pivi)' 



PoJ \ Pi 



(P2^2)' 



where 



Ao = l, A^ = {-irl[f,,. (5.17) 

A:=0 

\i y = Un {z) is a polynomial of degree n, then t>„ = yi"^ (2) = constant and we arrive at the 
Rodrigues formula for polynomial solutions of (15.21) . 

yn{z) = ^Acr''iz)piz)f''\ (5.18) 
P{z) 

where Bn = A~^yli''^ is a constant. These solutions correspond to the eigenvalues (15. lip . 

The polynomial solutions of (15. 2p obey an orthogonality property. Let us write equations for 
polynomials (x) and ym (x) in the self-adjoint form 

(cr (x) p (x) y'^ (x))' + A„p (x) yn (x) = 0, 

(a (x) p (x) (x))' + \mp {x) y^n {x) = 0, 
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multiply the first equation by (x) and the second by ?/„ (x) , subtract the second equality from 
the first one and then integrate the result over x on the interval (a, b) . Since 

Vm (x) {a (x) p (x) y'^ (x))' - y„ (x) {a (x) p (x) y'^ (x))' 
= ^ [cr {x) p {x) W {ym (x) , yn (x))] , 
where W {u, v) = uv' — vu' is the Wronskian, we get 

(Am - A„) / ym {x) Vn {x) p {x) dx = [a (x) p (x) W {ym (x) , (x))] |^^^ . (5.19) 

J a 

If the conditions 

a(x)p(x)x'=|^^^, = 0, A; = 0,1, 2, ... (5.20) 

are satisfied for some points a and b, then the right hand side of (15.191) vanishes because the 
Wronskian is a polynomial in x. Therefore, we arrive at the orthogonality property 

-6 

ym (x) yn (x) p (x) dx = (5.21) 

provided that 7^ A„. One can replace this condition by m 7^ n due to the relation A^, — A„ = 
(m - n) (r' + (n + m - 1) a"/2) if r' + (n + m - 1) a"/2 ^ 0. 

We shall refer to polynomial solutions of (15.21) obeying the orthogonality property (15.211) with 
respect to a positive weight function p (x) on a real interval (a, b) as classical orthogonal polynomials. 

Equation 

(a(x)p(x))' = r(x)p(x) (5.22) 

for the weight function p (x) is usually called the Pearson equation. By using the linear transfor- 
mations of independent variable x one can reduce solutions of (15.221) to the following canonical 
forms 

'l — x)°'{l + x)^ for cr(x) = l — x^, 



p[x) 



for a(x)=x, (5.23) 



„2 



e ^ for a (x) = 1. 



The corresponding orthogonal polynomials are the Jacobi polynomials P^' (x) , the Laguerre 
polynomials (x) and the Hermite polynomials Hn (x) . 

The basic information about the classical orthogonal polynomials is given in Table, which contains 
also the leading coefficients in the expansion yn{x) = a„x" + bnX^~^ + . . . , the squared norms 

dl= [ yl{x) p{x)dx (5.24) 



and the coefficients of the three-term recurrence relation 

X yn{x) = an yn+i{x) + (3^ yn{x) + 7„ yn-i{x) (5.25) 

with 

a. 



bn _ bn+1 _ d'^ 

) y-a ~ ~ 5 In ^n—l 



Pn= , In = -J^- (5.26) 
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Vnix) 


PlC'^Kx) (a >-!,/?>-!) 


Llix) (a > -1) 




{a,b) 
p{x) 

cr(x) 
r(x) 

Bn 

an 

K 

(p 

(Xn 

Pn 
'~1 n 


(-1,1) 

/?-«-(« + /? + 2) a; 
n{a + l3 + n + l) 

(-1) 
2"n! 

r(a + ^ + 2n + l) 

2'^ra!r(a + /5 + + 1) 
(a-/5)r(a + /3 + 2n) 


(0,oo) 

X 

1 + a — x 
n 

1 

(-1)" 


(— OO, oo) 

1 

-2x 
2n 

(-ir 

2" 


2"n!V? 

1 
2 




2"(n-l)!r(a + /3 + n + l) 
2«+/3+ir(a + n + l)r(/5 + n + 1) 


^ ' {n-\)\ 
r(a + n + 1) 


n!(a + /3 + 2n + l)r(a + /5 + n + 1) 

2(n + l)(a + /3 + n + l) 
(a + /3 + 2n + l)(a + /? + 2n + 2) 

(a + /3 + 2n)(a + /? + 2n + 2) 
2(a + n)(/3 + n) 

(a + /5 + 2n)(a + ^ + 2n + 1) 


n\ 

a + 2n + 1 
— (a + ra) 



More details about the Jacobi, Laguerre and Hermite polynomials and their numerous extensions 
can be found in [3j, ^42], [18], [55], [57], [TlJ, ^74j and references therein. 

5.3. Classical Orthogonal Polynomials and Eigenvalue Problems. The following theorem is 
a useful tool for finding of the square integrable solutions of basic problems in quantum mechanics 

Theorem 1. Let y = y{x) he a solution of the equation of hypergeometric type Ii5. 2\] and let p (x) , 
a solution of the Pearson equation Ii5. 22\) . he hounded on the interval (a, h) and satisfy the boundary 
conditions ^5.20\) . Then nontrivial solutions of ( 15. ^) such that y (x) a/ p (x) is hounded and of 
integrable square on (a, b) exist only for the eigenvalues given by Ii5.11\) : they are the corresponding 
classical orthogonal polynomials on (a, b) and can he found by the Rodrigues-type formula Ii5.18\) . 

The proof is given in [57j . 

5.4. Integral Representation for Special Functions. The differential equation of hypergeo- 
metric type (15.21) can be rewritten in self-adjoint form 

(a (z) p {z) y' iz))' + Xp (z) y (z) = 0, (5.27) 

where p (z) satisfies the first order equation 

{a{z)p{z)y = T{z)p{z). (5.28) 
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Nikiforov and Uvarov [57] suggested to construct particular solutions of the differential equation 
of hypergeometric type (15. 2p in a form of a general integral representation for special functions of 
hypergeometric type as a refinement of the Laplace method. A slightly modified version of their 
main theorem is 

Theorem 2. Let p (z) satisfy liS. 28\) and u be a root of the equation 

X + UT -l)a" = 0. (5.29) 

Then the differential equation ( (5.^) has a particular solution of the form 

f a'{s)p {s) 
p{z) J {s-z 



yi^) = yA^) = ^, l -r^^^ds, (5.30) 



c 

where is a constant and C is a contour in the complex s-plane, if: 

(1) the derivative of the integral 

v., (^) = / T^7%T ds wtth is) = (s) p (s) (5.31) 

can he evaluated for ji = u — \ and ji = u hy using the formula 

(2) the contour C is chosen so that the equality 



a^+\s)p{s) 



S2 



(5.33) 



[s-z) 

holds, where Si and $2 are end points of the contour C. 

We present here a simple proof of this theorem [71], which is different from one in |57j . 

Proof. The function (s) = (s) p (s) satisfy the equation 

(a(s)p,(s))' = r,(s)p,(s), (5.34) 

where Tj, (s) = r (s) + va' (s) . We multiply both sides of this equality by (s — zY''^^ and integrate 
over contour C. Upon integrating by parts we obtain 

(^{s)p^{s) 



{s-z) 



v+1 



S2 



Si 



(•'+l)rj^^''^= f'j^^"- (5.35) 

Ic [s - z) Jc [s - z) 



By hypothesis, the first term is equal to zero. We expand polynomials a (s) and (s) in powers of 
s — z : 

a (s) = aiz)+ a' {z) {s - z) + \a" (s - zf , 

Tu (s) = Ty {z) + rl {s- z) . 

Taking into account the integral formulas for the functions 99^ (z) , (2;) and (p^, ^_|_^ (z) , we 
arrive at the relation 

+ 1) U (^) u+l + ^' (Z) Vvv + Vv, u-l \ = {z) + < V'l., v-\- 
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Upon substituting t^, = r + ua' and using the formula y^'^^ = (z/ + 1) ^9^^ ^j^-^ one gets 

ay.;, + (a' - r) = (^r' + i (z/ - 1) a"^ y,,, (5.36) 

At the sane time, by differentiating the relation apy' = Cyaip^^^ we find that 

^apy' = ai^l^ + {a - r) <^^^. (5.37) 
Comparing (I5.36P and (I5.37P we obtain 

crpy' = K^C^ip^^ (5.38) 

where Kj, = r'+(z/ — f ) a" /2. Upon differentiating fl5.38p we arrive to the equation of hypergeometric 
type in the self-adjoint form 

{apy')' + Xpy = 0, 

where A = —pKu = —vt' — u {u — 1) a" jl. This proves the theorem. □ 

In the proof of Theorem 2 we have, en route, deduced the formula fl5.38p . which is a simple 
integral representation for the first derivative of the function of hypergeometric type: 

y'. {z) = I ds, (5.39) 



c 

(1) 



where Cl"' = k^C^ = ( r' + i (z/ - 1) a"j C^. Hence 



. / (5.40) 



where C^'^ = 17' Vr' + ^ (z. + p - 1) a"] C,. 

-l--t-p=o \ I J 



See [56], [7T] and [72] for an extension of this theorem to the case of the so-called difference 
equation of hypergeometric type on nonuniform lattices. 

5.5. Power Series Method. We can construct particular solutions of equation (15.21) by using the 
power series method; see, for example, the classical work of Boole [2T] . 

Theorem 3. Let a he a root of the equation a{z) = 0. Then Eq. / I5.^)) has a particular solution of 
the form 

oo 

y{z) = Y,(^n{z-ar, (5.41) 

where 

Cn+i ^ \ + n{T' + {n-l)a"2) 
Cn (n+ 1) (r(a) +na'(a)) ' ^' ' 

^f: 

d^ 

(1) lim — ym{x) = — y{x) with /c = 0, 1, 2; 
m-^oo ax dx 
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(2) lim (A - Xm) Cm{x - a)"" = 0. 



Em 1 
c„(x — a)" and Am = —mr' — -m{m — ■) 
n=0 2 

In the 
t{z) = 0, 



2 

In the case a{z) = constant ^ series { 5. 41 ) satisfies ( f5.^) when a is a root of the equation 



= __A±!^ (5.43) 
and convergence conditions (1 )-(2 ) are valid. 
Proof. The proof of the theorem relays on the identity 

where T^iO — '^iO + "^'^'(0 ^-^id A„ = — nr' — 72(7?, — l)o""/2, which can be easily verified. 
In fact, for a partial sum of the series (I5.41I) we can write 
/ d"^ d \ " 

( '^^^^J^ + ^^^^T ^ ^ I ^""^^^ ^ "^^^^ S ^"^^'^ - - «)""^ (5-45) 

^ ^ n=0 

m m 

+ ^ c„nr„_i (a) - a)"~^ + ^ c„ (A - A„) {z - a)". 

n=0 n=0 

By the hypothesis cr(a) = and the first term in the right hand side is equal to zero. Equating the 
coefficients in the next two terms with the aid of 

Cn+l A„ — A 



On (n + l)r„(a)' 

which is equivalent to f l5.42p . one gets 



(5.46) 



d^ d \ 

Taking the limit m ^ 00 we prove the first part of the theorem under the convergence conditions 
(l)-(2). 

When cr = constant we can obtain in the same manner 

(^{z)^ + r{z)^ + x]yU^) (5.48) 



dz"^ dz 

m m 

= cr ^ Cnn{n - l){z - a)""^ + ^ c„ (A - A„) {z - a)"" 

n=0 n=0 

= Cm (A — Am) {z — a)™, 

which proves the second part of the theorem in the limit m — 00. □ 
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Corollary. Equation Ii5. 2\] has polynomial solutions ym{x) corresponding to the eigenvalues A = 
Am = —mr' — m{m — l)a" /2, m = 0,1,2, .... 

This follows from flCTl) and fICTD . 

Examples. With the aid of linear transformations of the independent variable, equation (15.21) for 
r' 7^ can be reduced to one of the following canonical forms [SZ] : 

z{l - z)y" + [7 - (a + /? + l)z\ y' - a(3y = 0, (5.49) 

zy" + (7 - z)y' -ay = 0, (5.50) 
y" - 2zy' + 2vy = 0. (5.51) 
According to fl5.41l) - fl5.43l) the appropriate particular solutions are: 

the hypergeometric function, 

y{z) = ,F,{a, (3; 7; ^) = £ (5.52) 

^ (7)nn! 

the confluent hypergeometric function. 



(a)n 

yy^) = i-rH«,7, = 7 T 

n=0 ^ 

and the Hermite function. 



y{z) = ^F,{a; r,z)=Yl TVl ^^'^^^ 

n=o ^'y)^^- 



n=0 



2"r(i) ^ / 1 A 2'^r(-i) fi-u 3 2 



r(i^) V 2' 2' ; ' r(-|) V 2 '2' 

respectively. Here (a)„ = a(a + 1) . . . (a + n — 1) = r(a + n)/r{a) and r(a) is the gamma function 
of Euler. 

Generally speaking, these solutions arise under certain restrictions on the variable and parameters; 
see equation (18. ip below for more details. They can be extended to wider domains by analytic 
continuation. 

Extended Power Series Method. The solution (I5.41I) - (I5.42I) can be rewritten in the following 
explicit form 



where Cq is a constant. Using the expansion 



Cn+l ^a+n ~ A 

one can construct solutions of the more general form 

ra+fc(a)(a + k + 1) 



y(.) = 5^c„(.-0"+", ^ "I" , , (5.56) 

^ On (a + n + l)r«+„(a) 



n=0 A:=0 
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provided that a{a) = and aTa-i{a) = 0. In particular, putting a = we recover (I5.55p . 
We can also satisfy (15. 2p by using the series of the form 

if a{a) = and A = A-q,. Hence 

- aY ^ n (A - A_._._0 (. - a) 
When cr = constant 7^ one can write the solution as 

^ (X - a)°+'^ C„ A - \-a-n-2 

if r(a) = and A = A_q, (for even integer values of n) or A = A_a-i (for odd integer values of 
n). This method allows to construct the fundamental set of solutions of equation (15. 2p . Examples 
are given in [73j. Sec [9] and [72] for an extension of the power series method to the case of the 
difference equation of hypergeometric type on nonuniform lattices. 

5.6. Integrals for Hypergeometric and Bessel Functions. Using Theorem 2 one can obtain 
integral representations for all the most commonly used special functions of hypergeometric type, 
in particular, for the hypergeometric functions: 

2Fi(a,/3;7;^) = :pT-4r^ V T ^-^(1 - tr-"-i(l - zt)-^ cit, (5.61) 

r(a)r(7-a) Jo 

iFi(a;7;z) = — -^^^^ -l - ty-'^-^e'Ut, (5.62) 



r(a)r(7 — a) 



H,{z) = - — r / e-'-'^H-'-'dt. (5.63) 
[—^) Jo 

Here Re7 > Rea > and Re (— z^) > 0. These functions satisfy equations (I5.49p - (l5.5ip . respectively. 
The Bessel equation is 

z'^u" + zu + [z^ -u^)u = 0, (5.64) 

where z is a complex variable and parameter u can have any real or complex values. The solutions 
of (I5.64P are Bessel functions (z) of order u. With the aid of the change of the function u = (p{z)y 
when (p{z) = z^e^^ equation (15.641) can be reduced to the hypergeometric form 

zy" + {2iz + 2u+ l)y' + i{2u + 1)?/ = (5.65) 

and based on Theorem 2 one can obtain the Poisson integral representations for the Bessel function 
of the first kind, Jiy{z), and the Hankel functions of the first and second kind, Hj)^\z) and Hj)^\z) : 

H^y){z) = \ -— / e-*r-i/Ml±-) dt, (5.67) 

" ' \ nz r(z/+l/2) Jo \ 2zJ ' ^ ^ 
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where Rez/ > —1/2. It is then possible to deduce from these integral representations all the 
remaining properties of these functions. For details, see \57\ . 



The Lommel equation 



1 - 2a , ^ r/^-v'y-n2 - 



V H — V + 



v = 0, (5.68) 



where a, j3 and 7 some constants, is very convenient in applications. Its solutions are 

V (0 = Cu. (PC) , (5.69) 
where (z) is a Bessel function of order u. 

6. Solution of Dirac Wave Equation for Coulomb Potential 

This section is written for the benefits of the reader who is not an expert in relativistic quantum 
mechanics and quantum field theory. We separate the variables and construct exact solutions of the 
Dirac equation of in spherical coordinates for the Coulomb field. The corresponding four component 
(bispinor) wave functions are given explicitly by (I4.1ip - fl4.15p . We first construct the angular parts 
of these solutions in terms of the so-called spinor spherical harmonics or spherical spinors. 

6.1. The Spinor Spherical Harmonics. The vector addition j = 1 + s of the orbital 1 = —ir x V 

and the spin s = -cr angular momenta (in the units of h) for the electron in the central field gives 

the eigenfunctions of the total angular momentum j, or the spinor spherical harmonics |80j, in the 
form 

3^SW= E Ci-^^mJ^mAn) X.m. (6.1) 

mi+ms=m 

(j = |/-l/2|,/ + l/2; m = -j,-j + l, 
where Yim (n) with n = n{6,ip) = r/r are the spherical harmonics, Cj^ are the special Clebsch- 
Gordan coefficients, and ~ Xima^'^^ eigenfunctions of the spin 1/2 operator s : 

= ^X^m,, SsXim, = ^sXlm,, "1, = ±1/2 (6.2) 

given by 

' 







.0)' 


63], [55], and [80 


. From fferri) 






1/2 

/ J I, m- 
ms=-l/2 








Y 

I, m+ 




( — il 


1 ^1, r 
2 




1 /^im 


-2 ' 



X -I = U ; (6.3) 



Yl, m-ms (n) Xm, (6.4) 



.1 n 

2 



l=j± 1/2. 
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Substituting the special values of the Clebsch-Gordan coefficients [80j, we obtain the spinor spherical 



harmonics (n) = y,-^"*^^^^ (n) in the form 



/ _ /(j + l/2)T(m-l/2) 



(n) 



2j + (1 ± 1) 



^■±1/2, m-1/2 (n) 



'(j + l/2)±(m + l/2) 



V 



2j + (1 ± 1) 



^■±1/2, m+1/2 (n) 



(6.5) 



with the total angular momentum j = 1/2, 3/2, 5/2, ... and its projection m = —j, — j + j — 1, j. 



The orthogonality property for the spinor spherical harmonics (n) = 3^]:^^^^'* (n) is 



3^],^, (n) dw = Sii>6ii>6. 



(6.6) 



with duj = sin 6 dOdip and O<0<7r,O<(/9< 27r. They are common eigenfunctions of the following 
set of commuting operators 



fyfm (n) 



1 



1 



y%, (n) = J (j + 1) yf^ (n) 



hyfn, (n) 



"^3^im (n) 



1 3^ (n) 



j±- j±- + l 3;,t(n) 



But 



or 



33^f^ (n) . 



(6.7) 
(6.8) 
(6.9) 
(6.10) 



1^ + ^-1 + 



(6.11) 



This implies that the spinor spherical harmonics yf^ (n) are also eigenfunctions of the operator 



o-l: 



and it is a custom to write 



{a- ■ 1) yf^ (n) 



1 ± ( J + 2 ) ) ^ 



(n) , 



^■I)3^f„(n) = -(l + «:)3^; 



n 



(6.12) 



(6.13) 



1 



whe.e the quantum number . = = ± + - j = ±1, ±2. ±3, ... takes all positive aud negative 
integer values with exception of zero: k ^ 0. 

Finally, the following relation for the spinor spherical harmonics, 

(^ ■ n) 3^f^ (n) = -3^^ (n) , (6.14) 

plays an important role in the Dirac theory of relativistic electron. In view of (cr ■ n) = 1, it is 
sufficient to prove only one of these relations, say 

{a . n) y;^ (n) = -yr^ (n) , 
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and the second will follow. A direct proof can be given by using the recurrence relations for the 
spherical harmonics fl8.17p - fl8.19p . or with the help of the Wigner-Eckart theorem; see [HS] and 
the reader can work out the details. 



The quadratic forms 



t 



Qjm=[yfl{n)) y%{r,) (6.15) 



of the spinor spherical harmonics yf^ (n) describe the angular distributions of the electron in states 



■ jm 



with the total angular momentum j, its projection m and the orbital angular momentum /. These 
forms, given by [80] 



1 / I |2 I |2\ 

Qjm {Q) = ip ({j + m) |y,_i/2, m-1/2 (n) | + (j - m) |Fj_i/2, m+1/2 (n) | j (6.16) 



2j 

1 / I |2 I 

= 2j + 2 V"'' + + ^) l^i+1/2, m+1/2 (n) I + (j - m + 1) \ Yj+i/2, m-i/2 (n) 

are, in fact, independent of I and ^p. There is the useful expansion in terms of the Laguerre polyno- 
mials 

Q,m {0) = 5^ a. (j, m) P2s {cosO) (6.17) 

with the coefficients of the form 

M^-'-^^^m^) { ^^^^^ ^^^^^ } c;:vjo%oq:2.o 



- r_i 1^ / (2J + 2S + 1) (2j-2s)! (j + .-l/2)!(2s)! 

^ 47r V (2j + l)(2j + 2.)! (j - s - 1/2)! (s!)^ ^ ^ 

See [80j for more information. 

6.2. Separation of Variables in Spherical Coordinates. Using the explicit form of the a and j3 
matrices (14.30 we rewrite the stationary Dirac equation (14.90 in a central field with the Hamiltonian 

= + »c=^ + U(r)=(" ,j ) (6.19) 



and the bispinor wave function 



^ = ( ^ ) (6.20) 



in a matrix form 



or 



^ + 2 V '')=e(^ (6.21) 

ccrp t/ — mc J \ X J \ X 



cap if = [E + mc^ - U) x, (6-22) 

cap X = [E-mc^ - U) ip. (6.23) 
Here we shall use the following operator identity 

cr- V = (cr-n)(n- V + io--(nxV)) (6.24) 



36 SERGEI K. SUSLOV AND BENJAMIN TREY 

in the form 

cap = he (crn) f -nV + -crl ) , (6.25) 
\z r J 

where 1 = —ir x V is the operator of orbital angular momentum, n = r/r and p = —ih^. It can 



be obtained as a consequence of a more general operator identity |31j| 

(cr ■ A) (cr ■ B) = A • B + i<T- (A X B) , (6.26) 

which is valid for any vector operators A and B commuting with the Pauli cr-matrices; it is not 
required that A and B commute. The proof uses a familiar property of the Pauli matrices 

a, o-fc = iciki oi + 5ik, (6.27) 

where Ciki is the completely antisymmetric Levi-Civita symbol, 5ik is the symmetric Kronecker delta 
symbol and we use Einstein's summation rule over the repeating indices; it is understood that a 
summation is to be taken over the three values of / = 1, 2, 3. Thus 

(cr ■ A) (cr ■ B) = (a, A,) (fXfc B^) 
= {ai CTfc) Ai Bk = iai Ai Bk + 5ik Ai Bk 
= iai (A X B), + Ak Bk = icr- (A x B) + A ■ B, 

where (A x B)^ = euk Ai Bk = Ciki Ai Bk in view of antisymmetry of the Levi-Civita symbol: 

^iki = —e-iik = Giik, and 6ik Ai = Ak- 

If A = B, Eq. (16.261) implies (cr ■ A)^ = A^. In particular, (crn)^ = = 1, and the proof of the 
"gradient" formula fl6.24p is 

cr ■ V = (cr ■ n)^ (cr ■ V) = (o- ■ n) ((cr • n) (cr ■ V)) 
= (cr • n) (n • V + icr ■ (nx V)) 

by (16.261) with A = n and B = V. In a similar fashion, one can derive the following "anticommu- 
tation" relation, 

(crn) (crl) + (crl) (crn) = -2 (crn) , n = r/r, (6.28) 
we leave details to the reader. 

The structure of operator crp in f l6.25p suggests to look for solutions of the Dirac system fl6.22p - 
(I6.23P in spherical coordinates r = r n{6,(p) in the form of the Ansatz: 

cp = ip(r)=yin) F(r), (6.29) 
X = x{r) = ~^ ((^rn) y (n)) G (r) , (6.30) 

where y = y^{n) are the spinor spherical harmonics given by (16.50 . This substitution preserves 
the symmetry properties of the wave functions under inversion r — > — r. Then the radial functions 
F (r) and G (r) satisfy the system of two first order ordinary differential equations 

^ + i±^F = + G, (6.31) 

dr r he 

^ '^-'-E + Uir)^^ (6.32) 

dr r he 
where k, = k± = ± {j + 1/2) = ±1, ±2, ±3, ... , respectively. 
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If/ = /(r)=/(rn),then 

a/ ^ d£dr ^ 
dr dr dr 

and in spherical coordinates Eq. f l6.25p becomes 

ccrp = he (crn) | + -crl ) . (6.33) 

\i or r J 



Thus 



fid i \ 
ccrp (p = he (crn) -— + -cr\ yP 

\t or r J 



= -ihe{(Tny) — + ^—F 

\dr r 

by (16.131) . and we arrive at (16.311) in view of (16.221) and (16.301) . Equation (16.321) can be verified in a 
similar fashion with the help of (16.281) or (I6.14p . 

Eqs. (I6.31I) - (I6.32I) hold in any central field with the potential energy U = U {r). For states with 

discrete spectra the radial functions rF (r) and rG (r) should be bounded as r — > and satisty the 
normalization condition 

J^^ ^ dv = J r^ {F^ (r) + (r)) dr = 1 (6.34) 
in view of (gTD, dOOjl - dCT]) and M . 



6.3. Solution of Radial Equations. For the relativistic Coulomb problem U = —Ze'^/r, we 

introduce the dimensionless quantities 

E me Ze^ 

X = pr = ^— r, fi = — — (6.35) 



mc^' h ' he 

and the radial functions 

/(x) = F(r), g{x) = G{r). (6.36) 
The system (I6.3ip -( l6.32p becomes 

df + K 

dx X '' \ ' ' X 

= (l-,-t);^ (6.38) 



/ = [l + E+f^)g, (6.37) 



dx X \ X/ 

We shall see later that in nonrelativistic limit c ^ oo the following estimate holds |/ (x)| ^ \g {x) \ . 

We follow with somewhat different details. Let us rewrite the system (I6.37p - (l6.38p in matrix 
form [13]. If 

^ ^ I "1 \ ^/ (^) ^ u' =1'''^ \ (6 39) 

U2 ] \ xg{x) j ' \ 
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Then 

u' = Au, (6.40) 

where 





( an 


ai2 ^ 












^ 021 


^22 y 






(6.41) 

To find u\ [x) , we ehminate U2 [x] from the system fl6.40p . obtaining a second order difi^erential 
equation 

u[ - f an + a22 + — ) A (6-42) 

V ^12 / 

( I ^12 \ 

+ 011^22 — 0,120.21 — On H On ] Ui = 0. 

V Oi2 / 

Similarly, eliminating Ui [x) , one gets equation for U2 (x) : 



a'r, 



ul - On + a22 + — «2 (6.43) 

021, 



/ / O21 \ 
+ 011022 — 012O21 — O22 H 022 ^2 = 0. 

V «21 / 



The components of matrix A have the form 

aik = hik + Cik/x, (6.44) 

where hik and are constants. Equations fl6.42p and fl6.43p are not generalized equations of 
hypergeometric type (15. ID . Indeed, 



o'i2 C12 



012 Ci2a; + 6i2x2' 
and the coefficients of u'^ {x) and Ui {x) in (16.420 are 

oi2 Pi {x) C12 



On + 022 + 



O12 X C12X + hi2x'^ ' 

, a'i2 P2 {x) C12 (cn + hiix) 



011O22 — 012O21 — a^^ H On — — 5 — — 5-, 

012 x'^ (ci2 + 012a;) X'^ 

where pi (x) and p2 (x) are polynomials of degrees at most one and two, respectively. Equation 
(I6.42P will become a generalized equation of hypergeometric type (15. ip with a (x) = x if either 
612 = or C12 = 0. The following consideration helps. By a linear transformation 

V2 j \ "2 y 

with a nonsingular matrix C that is independent of x we transform the original system (16.40p to a 
similar one ^ 

v' = Av, (6.46) 

where 

vi \ ~ 1 / On ai2 

A = CAC-^ 
V2 I \ 021 O22 
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The new coefficients are linear combinations of the original ones aj^. Hence they have a similar 
form 

= bik + Cik/x, (6.47) 

where bik and Cik are constants. 

The equations for vi {x) and V2 {x) are similar to fl6.42p and (I6.43p : 

v'l - (au + + ^ ) v[ (6.48) 

+ 011022 - 012021 - Oil + — an Vi = 0, 
V ai2 / 

t'a - ( «ii + ^22 + ^^2 (6-49) 

V «2l/ 



+ 011022 - «12a21 - 022 + — '^22 "^2 = 0. 
V «21 / 



The calculation of the coefficients in (16.481) and (16.491) is facilitated by a similarity of the matrices 
A and A : 

an + 022 = 0-11 + O22, 0-110.22 ~ O12O21 = O11O22 "~ O12O21. 

By a previous consideration, in order for ( 16.481) to be an equation of hypergeometric type, it is 
sufficient to choose either 612 = or C12 = 0. Similarly, for ( 16. 49^ : either 621 = or C21 = 0. These 
conditions impose certain restrictions on our choice of the transformation matrix C. Let 

C=f!'n. (6.50) 



Then 



and 



7 S 



C-^ = ^ ^ ) , A = det C = a5 - /37, 



A \ —7 a 



A = CAC^^ (6.51) 

1 / oiia(5 — 0120:7 + 021/95 — 022/?7 0120^ — 012/^^ + (022 — ou) «/5 

A \^ a2i5^ - Oi27^ + (oii - O22) 7*^ 01207 - oii/?7 + a22aS - a2i(35 
For the Dirac system (lOnil - dOTD : 



ail = , O12 = 1 + e + -, 

X X 

_ 1 ^^ _ 
O21 — 1 — 6 , O22 — — 

X X 



and 



A ai2 = a^-P^+ {a^ + P^) e + ^ — , (6.52) 

A ~ X2 2 /r2 , 2^ {S^ + -f^) fi + 2-f6K 

A 021 = 5 - 7 - [5 + 7 j e - ^ . (6.53) 
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The condition 612 = yields (1 + e) - (1 - e) /J^ = 0, 
£12 = " (a^ + + 2apK = 0, 
621 = " (1 + 5)72 -(1-6)52 = 0, 
C2i = " + 72) /i + 27(5fi: = 0. 

We see that there are several possibilities to choose the elements a, (3, 7, S of the transition matrix 

C. All quantum mechanics textbooks use the original one, namely, 612 = and 621 = 0, due to 
Darwin ^30] and Gordon |13]; cf. equations fl6.83p - fl6.84p below. Nikiforov and Uvarov [57j take 
another path, they choose C12 = and C21 = and show that it is more convenient for taking the 
nonrelativistic limit c 00. These conditions are satisfied if 



where u = a/k^ — , and we finally arrive at the following system of the first order equations for 
f 1 (x) and f 2 (x) : 

< = (—- z) ^1 + f 1 + — ) V2, (6.55) 















) + 


(^ 


V 










Vx + 1 




1 - 


V ) 


(- 



V 



v'2 = {l~—)vi+(--^-^\v2. (6.56) 



Here 



Tr A = an + 022 = 0, det A = e'^ - I + — - z/^ = _ (6.57) 

X x^ 

which is simpler that the original choice in [ST]. The corresponding second order differential equa- 
tions ([S3HD-(IS3nD become 

x^ 

,,^^ ^-'-^)-'^^-^---^--^) = 0. (6.59) 
x^ 

They are the generalized equations of hypergeometric type (15. ip of a simplest form r = 0, thus 
resembling the one dimensional Schrodinder equation; the second equation can be obtain from the 
first one by replacing v ^ —v. 

Let 1 + EKjv = 0, then e = —u/k that is possible only if /t < 0, since u > and e > 0. The 
corresponding solution of (16.551) . 

vi (x) =ClX-'^e("'^^)/^ 
satisfies the conditions of the problem only if Ci = 0. Then from fl6.56p 

(x) =C2X^e-("'^^)/^ 

which does satisfy the condition of the problem with C2 7^ 0. 

Let us analyze the behavior of the solutions of fl6.58p as x ^ 0. Since 

I (e^ _ 1) + 2£/ix| < 1/ (z/ + 1) 

as X — > 0, one can approximate this equation in the neighborhood of x = by the corresponding 
Euler equation 

x'^v" — z/ (z/ + 1) t>i = 0, 



THE HAHN POLYNOMIALS IN THE COULOMB PROBLEMS 



41 



whose solutions are 

vi (x) = Cix^+^ + C2X-^ C2 = 0. 

Thus vi — > Cix^"*"^ as X — 0. The resuhs for f l6.59p are similar: V2 C2X" as x ^ 0; one can use 
the symmetry v ^ —v. 

Equation (16.581) is the generalized equation of hypergeometric type (15.11) with 

a (x) = X, r (x) = 0, 

a (x) = [s^ - 1) x^ + 2enx - z/ (z/ + 1) . 



The substitution 



where 



vi=^{x)y{x), — = —-. 

(f a [X) 



a — T 



71 X 



± 



a — r 



a 



kcr 



(6.60) 
(6.61) 



with k = \ — it' and r (x) = r (x) + 27i (x) , results in the equation of hypergeometric type 

a (x) y" + T{x)y' + Xy = (6.62) 
by the method of [57|; see also Section 5.1. From the four possible forms of tt (x) : 



71 (x 



) = i±(v/r3Fx±(.+i)) 



(6.63) 



corresponding to the values of k determined by the condition of the zero discriminant of the quadratic 
polynomial under the square root sign in (16.611) : 

k-2eii = ±Vl - £2 {2u + 1) , (6.64) 

we select the one when the function r (x) has a negative derivative and a zero on (0, +00) . This is 
true if one chooses 

k = 2ei2-a {2v + 1) , 
TT (x) = V + 1 — ax, 
r (x) = 27r (x) = 2 (z/ + 1 — ax) , 

A = k + 7i' = 2{e^- a{v +1)) 



and 

where a = \/l — and u 
u ^ -u in ([6l3])-([6JlD. 



^ (x) = x"+^e-"", p (x) = x'^+^e""'^", 



a/k^ — p'^. The analysis for (16.591) is similar, one can use the symmetry 



From (ESI and (109l) 

{F^ (r) + (r)) dr = (3^^ 



{u\ (x) + u\ (x)) dx = 1. 



(6.65) 



It requires by (I6.45P the square integrability of vi (x) and V2 (x) . Their boundness at x = follows 
from the asymptotic behavior 0. So 



vl (x) dx 



(x) (x) (ix 



X y'^ (x) p (x) dx < 00. 



(6.66) 
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For the time being, we replace this condition by 

r>00 

(x) p (x) dx < CO (6.67) 



in order to apply Theorem 1, and will verify the normalization condition fl6.65p later. Then the 
corresponding energy levels e = En are determined by 

A + nr' + (n - 1) a" = (n = 0, 1, 2, ... ) , (6.68) 

whence 

e/i = a (z/ + n + 1) , (6.69) 
and the eigenf unctions are given by the Rodrigues formula 

Vn (x) = 4\ (^) P = Cn x-^'^-'e^'^-f^ [x'^^-^'e-^'^^) . (6.70) 

p [x) ax"- 

The functions ?/„ (x) are, up to certain constants, the Laguerre polynomials L'^^~^^ (0 with ^ = 2ax. 

The previously found eigenvalue e = —u/k satisfies f l6.69p with n = —1. Consequently it is natural 
to replace n by n — 1 in fl6.69l) - fl6.70p and define the eigenvalues by 

6p = a{u + n), a = Vl-e^ (n = 0, 1, 2, ... ) . (6.71) 

Solving for e gives the Sommerfeld-Dirac formula fl4.16p . The corresponding eigenfunctions have 
the form 

f 0, n = 0, 

^''<^>n^.r«e-^/^eV(f). 1.2.3..... ''"'^ 

They are square integrable functions on (0, oo) . The counterparts are 

V2 (x) = BXe-^^'Ll^-' (0 , n = 0, 1, 2, ... . (6.73) 

It is easily seen that our previous solution for e = —u/k, is included in this formula when n = 0. By 
equation (16.551) the other solutions can be obtain as 

{x) = ^^-i^ {x) + - f ) ^1 (x)) 

and substituting vi (x) from (I6.72p one gets 

V2 (x) = Ce~^"Y (0 , 

where Y {S) is a polynomial of degree n. But function f 2 (x) satisfies fl6.59p . By the previous 
consideration the substitution 

gives 

xy" + (2z/ — 2ax) y' + 2any = 0, 

in view of the quantization rule (16.711) . The change of the variable y (x) = Y (C,) with ^ = 2ax 
results in 

^Y" + (2z/ -OY' + nY = (6.74) 

and the only polynomial solutions are the Laguerre polynomials L"^'^ (0 , whence (16.731) is correct. 
Solutions f 2 (x) are square integrable functions on (0, 00) . 



V2 (x) = x^e "^2/ (x) 
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To find the relations between the coefficients An and i?„ in fl6.72p and fl6.73p we take the hmit 
X — > in (16.551) with the help of the following properties of the Laguerre polynomials [SJ], [S7|, |74j : 

^l:, (0 = io , (0) = ^^,^+"+1^ (6.75) 

at, nlT [a + 1) 

The result is 

2a {u + 1) A„LtV (0) = -2auAnLl^_\' (0) + (l + ^) i?nC~' (0) , 

whence 



Since 



An= "".^^1 . {n = 1,2,3, ...). 

an [n + 21/) 



a^n{n + 2u) = a^ [{n + vf - v^) = - v'^ 



2\ 2 2 2 2 



we have proved the useful identity 

a^n{n + 2v)=e^K^ -v^, (6.76) 

and the final relation is 



a 



An = Bn. (6.77) 

— V 



By ([S3SD and (1^31) we find 



Therefore 



M2 y \ ^2 / 2Z/ (k - I/) y K - 

^/(^) = ^-Tf^re-^/^M^n-V (0 + Mr MO) , (6.78) 
X9{x) = ^^Y^re-^/^^i^C-V (0 + ^72Lr MO) , (6.79) 

. o/^ /. a{K-iy) 
fi = , f2 = t^-v, gi = , g2 = /i. (6.80) 

SK — V EK — U 

These formulas remain valid for n = 0; in this case the terms containing L'^_^i (0 have to be taken 
to be zero. Thus we derive the representation for the radial functions (I4.13P up to the constant Bn- 
The normalization condition (16.650 gives the value of this constant as 



where 



5. = a/?^/%/i^l^4^^^1^. (6.81) 

This has been already verified in Section 4.5. Observe that Eq. (16.811) applies when n = 0. 

The familiar recurrence relations for the Laguerre polynomials (I8.5p - (I8.6I) allow to present the 
radial functions (14.131) in a traditional form [1], [16], [28] as 



Fir) 



G{r) ] V ij{K-iy) {en -u)T{n + 2v) 
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X 



Pi 



a2 
P2 



iO 



with 



(6.82) 



(6.83) 



0^2 = — vT+T ^(k — z/) Vi + £ — /iVi — £^ 

/?2 = Vl - £ ((/t - z^) Vl + e - /iVl - . (6.84) 
By (18.21) one can rewrite this representation in terms of the confluent hypergeometric functions. 



6.4. Nonrelativistic Limit of the Wave Functions. Throughout the paper we have always 
used the notation n = rij. for the radial quantum number, which determines the number of zeros of 
the radial functions in the relativistic Coulomb problem; see f l4.13p . For the sake of passing to the 
limit c — > cxD in this section, let us introduce the principal quantum number of the nonrelativistic 
hydrogen atom as n = + |k| = nr+j + 1/2 and temporarily consider = n^ + z/ as its "relativistic 
analog" . As c — > 00 one gets 



V 



2 -/i2 



2\k\ 8 



\K 



2|/t| 8|k|' 

as /i = Z(? jhc — > 0. As a result, for the discrete energy levels 

we arrive at the expansion (14.171) in the nonrelativistic limit c 
In a similar fashion, 



^ +0(/.^), 



00. 



1 + 



/i 



as — >• 0, thus giving 
by (EJl). Also 



mc 
2a — — r 
n 



2 |k| (rij. + |k|) 



+ O (/) 



(i + o(/.^)) 



2Z f r 

lim qyc) = rj = — — 



2 k 



£K — Z/ = (k — |k|) + 



(^r 


+ 


I'^l)'- 


- K, 




2 




(n^ + 





as /i ^ 0. This allows to evaluate the nonrelativistic limit of the transition matrix: 



S 



fi 
91 



/2 
92 



f 



a/i 



EK — V 

a{K — u) 

£K — U 



K — V 



\ 



(6.85) 
(6.86) 

(6.87) 



(6.88) 
(6.89) 

(6.90) 

(6.91) 
(6.92) 

(6.93) 
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There are two distinct cases with the end result 

y^F \ ( ±y^R 
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± 







/i ^ 0. 



(6.94) 



Here R = Rni (r) are the nonrelativistic radial functions fl3.2p - fl3.3l) and 3^^ = 3^]^^'^^'' (n) are the 
spinor spherical harmonics (16. 5p . 



Indeed, if k = |k| = j + 1/2 = /, 

/ 2k, [rir + K 

S = S+ (fi) = 



rir [rir + 2k 
{rir + k) jj, 

\ rir {fir + 2k) 



+ 0(^2 



2 k 



+ 0(/i3 



+ 0(/i^ 



/i /i 



as /i or 



1 



lim 5"+ (/i) = ,„ . 



In this case u — » / and, therefore. 



F(r) 
G(r) 



Ze^m 



3/2 



(n-Z- 1) 
(n + Z)! 



in the limit c oo thus giving 







/i 



0. 



(6.95) 



(6.96) 



(6.97) 



In a similar fashion, when k 



\K\ 



- (j + 1/2) = — / — 1 one gets 

ly^G y V 

due to the corresponding asymptotic form of the transition matrix S = (/i) : 



/i ^ 



( 



\ 



2 |k| (rir + |fi:|) 



+ 0(/i4) -2|ft:| + 0(/i2) \ 



0(/x3) 



as /X — >^ [57]. This completes the proof of (I6.94p . 



1 



(6.98) 



(6.99) 



The representation of the radial functions in the form fl4.13p . due to Nikiforov and Uvarov |57j . 
is well adapted for passing to the nonrelativistic limit since one coefficient of the transition matrix 
S is much larger than the others as /x — 0. In the traditional form f l6.83p -( l6.84p . however, there is 
an overlap of the orders of these coefficients and one has to use the recurrence relations f l8.5p -( l8l6l) 
in order to obtain the nonrelativistic wave functions as a limiting case of relativistic ones. 
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7. Method of Separation of Variables and Its Extension 

In this section we give an extension of the method of separation of variables, that is used in 
theoretical and mathematical physics for solving partial differential equations, from a single equation 
to a system of partial differential equations which we call Dirac-type system. 

7.1. Method of Separation of Variables. We follow [57] and give an extension for suitable 
Dirac's systems. The method of separation of the variables helps to find particular solutions of 
equation 

£m = (7.1) 
if the operator L can be represented in the form 

L = MxMx^ M2M2. (7.2) 

Here the operators M.\ and M.2 act only on one subset of the variables, and the operators M\ and 
M2 act on the others; a product of operators M-iMk means the result of applying them successively 
{M.iNk)u = J^i{Afku) with 2, /c = 1, 2 ; it is assumed that the operators Aii and A/i are linear 
operators. 

We look for solutions of equation (17.11) in the form 

u = fg, (7.3) 

where the first unknown function / depends only on the first set of variables and the second function 
g depends on the others. Since 

MMu = {M-M{f g) = Md^k{f g)) 
= if {Mug)) = (MJ) {Mug) 
the equation Cu = can be rewritten in the form 

Mif ^ Af2g 
M2f Afig' 

where the left hand side is independent of the second group of the variables and the right hand side 
is independent of the first ones. Thus, we must have 

Mif ^ _U2g_ ^ 
M2f M^g ' 
where A is a constant, and one obtains equations 

Mif = XM2f, U2g = -\Nig (7.4) 

each containing functions of only some of the variables. Since £, is linear, a linear combination of 
solutions, 

u = ^Ckfkgk (7.5) 

k 

with some constants Ck, corresponding to all admissible values of A = Afc, will be a solution of 
the original equation (17.11) . Under certain condition of the completeness of the constructed set of 
particular solutions, every solution of (17.11) can be represented in the form (17. 5p . The method of 
separation of variables is very useful in theoretical and mathematical physics and partial differential 
equations — including solutions of the nonrelativistic Schrodinger equation — but, as we have seen 
in Section 6.2, it should be modified in the case of the Dirac equation. 
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Example. The Schrodinger equation in the central field with the potential energy U (r) is 

2^77 

A^ + — {E-Uir))ij = 0. (7.6) 
The Laplace operator in the spherical coordinates r, 6, ip has the form [55], [57] 

A = A, + ^A^ (7.7) 



with 



Thus 



2m 1 
A^i = A, + — (i?-f/(r)), M2 = ^ (7.9) 

M=id = /, Af2 = A^ (7.10) 
and separation of the variables = R{r)Y {6, ip) gives 

A^Y{e,v) + XY{e,^) = 0, (7.11) 

Bounded single- valued solutions of equation (17.111) on the sphere S"^ exist only when A = / (/ + 1) 
with I = 0,1,2, .... They are the spherical harmonics Y = Yi^ [9, ip) . 

7.2. Dirac-Type Systems. Let us consider the system of two equations 

Vu = av, (7.13) 

Vv = (3u, (7.14) 

where u = u (x) and v = v (x) are some unknown (complex) vector valued functions on R" (or 
) . Here operator V has the following structure 

V = U{n){Vi{r)Ci{n)+V2{r)C2{rv)), A^^ (n) = id = /, (7.15) 

where Vi = Vi (r) , Ck = C-k (n) and N = N (n) are linear operators acting with respect to two 
different subsets of variables, say "radial" r and "angular" n variables, respectively (in the case 
of the hyperspherical coordinates in R"' [55] one gets x = rn and = 1, which justifies our 
terminology). The following algebraic properties hold 

[V„Ck\ = [Vi,M]=0, (7.16) 

[M,C^] = [Cr,C2] = 0, (7.17) 

MC2 + C2M = lAf, (7.18) 
where [A, B] = AB — BA is the commutator and 7 is some constant. 
We look for solutions of fl7.13l) - fl7.14l) in the form 

u = y{n)R{r), (7.19) 

V = i^^y (n)) s (r) , (7.20) 

where y is the common eigenfunction of commuting operators Ci and C2 '■ 

Ciy = K,y, C2y = K2y. (7.21) 
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If w = w (x) = F (n) G (r) , we define the action of tfie "radial" and "angular" operators in (17.151) 
as follows 

dw = (CiF) G, Mw = (AfF) G, VkW = F (VkG) . (7.22) 
The Ansatz (I7.19l) -( 17.20l) results in two equation for our "radial" functions R and S : 

K1V1R + K2V2R = aS, (7.23) 
K^V^S + {-f-K2)V2S = PR. (7.24) 

Indeed, in view of f l7A3l) -f l7l4l) and (17:191) -( 17:211) one gets 

Vu = U{ViCi+V2C2)yR 

= U{{Ciy) {v,R) + {C2y) {V2R)) 

= {My){K{D^R + K2V2R) 

= a{^y)S = av, 
which gives (I7.23p . In a similar fashion, with the aid of (17.181) 

Vv = U{ViCi + V2C2){Uy)s 

= Af {{CiAfy) (v^s) + {C2My) {V2S)) 

= {MCiUy) (ViS) + {^C2^y) {V2S) 

= {M'c,y) {v,s) + ((7 - C2) M^y) {V2S) 

= y{K{D^s + {n-K2)V2S)=(3yR = i3u, 

which results in the second equation (I7.24p and our proof is complete. 
Example. The original Dirac system (I6.22I) - (I6.23I) has 

^ U \ [if 

and 



Here 



and 



V = ccT^ = tu^{a^)[\^^ + '-a\]. (7.25) 

fic d 1 f%c 

Af=cTn, Vi = ——, £i = id = /, V2 = — , £2 = 0-1 (7.26) 

t or r 



a{r) =E + mc^ -U (r) , (3 (r) = E - mc^ - U (r) , 7 = -2 (7.27) 

by (16.281) . Moreover, ki = 1, K2 = ^ (1 + ^) and we use R = F (r) , S = —iG (r) . The Anzats 
(I6:29|) - (l6:30|) gives the familiar radial equations (lOTD - (16321) . 

8. Appendix: Useful Formulas 

This section contains some relations involving the generalized hypergeometric series, the Laguerre 
and Hahn polynomials, the spherical harmonics and Clebsch-Gordan coefficients, which are used 
throughout the paper. 

The generalized hypergeometric series is [3], [IDj, [22], [12] 

pFq{ai, a2, ... , ap] bi, 62, ••• , bq] z) (8.1) 
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p I S \ _f.(ai)„(a2)„...(a,)„ 



^1' ^ (&l)J&2)„...(Mn^! 



where (a)^ = a (a + 1) ... (a + n — 1) = F (a + n) /F (a) . By the ratio test, the pFq series converges 
absolutely for all complex values of z if p < g, and for \z\ < 1 if p = g + 1. By an extension of the 
ratio test ([21], P- 241), it converges absolutely for \z\ = 1 ii p = q + 1 and z^()oip = q + l and 
Re [61 + . . . + 6q — (ai + . . . + Op)] > 0. If p > g + 1 and z^Qoip = q + l and \z\ > 1, then this 
series diverges, unless it terminates. 

The Laguerre polynomials are defined as [3], [55], [57], [71] 



F(a + n + 11 



-72 



It is a consequence of Theorem 3. The differentiation formulas [55], [57] 

^L^J (x) = -L^^ll (x) , (8.3) 

(x) = nL- (x) - (a + n) L^^ (x) (8.4) 

imply a recurrence relation 

xLltl (x) = (« + n) (x) - (x) . (8.5) 
The simplest case of the connecting relation fl2.12p is 

LI (x) = L^^ (x) - L^il (x) . (8.6) 

The Hahn polynomials are [55], [57] 

^" ^^^^)-^-^) n\T{N-n) /? + l, 1 - iV ' " ^^'^^ 

We usually omit the argument of the hypergeometric series 3F2 if it is equal to one. An asymptotic 
relation with the Jacobi polynomials is 

^ (j (1 + ^) - = is) + O (^) , (8.8) 

where N = N + {a + (3) /2 and N ^ 00; see [55] for more details. 
Thomae's transformation [101, US) is 



-n, a, 6 \ _ h\ ( —n, c — a, b 



c, d ' " / (d)^ ^"M c, 6-d-n + l 



3^2 I ^ , i 1 I = 3F2 I, , . ^ 1 ; 1 (8-9) 



with = 0, 1,2, ... . 

The summation formula of Gauss [3], [10], [12] 



, ^ \ F(c)F(c-a-6) , , , 

2F1 ^ ; 1 = TV, Re(c-a-6)>0. (8.10) 



c ' / T{c-a)T{c-b)' 
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The gamma function is defined as [3], [36], [57] 

/>oo 

T{z)= / e"V~^ dt, Re z>0. (8.11) 



It can be continued analytically over the whole complex plane except the points z = 0, —1, —2, ... 
at which it has simple poles. Functional equations are 

T{z + l)=zT{z), (8.12) 
T{z)T{l-z) = -^, (8.13) 

Sm TTZ 

2^'-'T {z)T{z + 1/2) = V^T{2z). (8.14) 

The generating function for the Legendre polynomials and the addition theorem for spherical 
harmonics give rise to the following expansion formula [57], [80] 

1 °° ' 47r 

^t-th = E E ^nfi ^™ (^1' v^i) (02, V,) , (8.15) 

where r< = min (ri, r2) and r> = max (ri, r2) . 

The Clebsch-Gordan series for the spherical harmonics is [55], [63], [80] 

Yi,mA0,^) Yhm2{0,^) = 2^ W , . X (8.16) 

where C\^^i^j^^ are the Clebsch-Gordan coefficients. The special case /2 = 1 reads [39] 



-sin^e- Y, = / (^ + y,^, ^ - J ' ^^7 T ~ !^ ^'-i (8-17) 
' ^' (2/+l)(2/ + 3) ^' V (2/ + l)(2/-l) V 



■n-^^v (/ - m) (/ - m + 1) / (/ + m)(/ + m + l) 

sm^e n ..+1 = W (2/ + i)(2/ + 3) ™ ~ V (2/ + 1) (2/ - 1) ^^'^^^ 



where 



= ^/(2/+l)(2/ + 3) - + V (2/ - 1) (2/ + 1) 
87r , • I Ait 

— Fi, ±1 = Tsin^e^^^ V Y = ^^-^^^ 
These relations allow to prove 06.141) by a direct calculation. 
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